Fractional Excitations in the Luttinger Liquid. 
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We reconsider the spectrum of the Luttinger liquid (LL) usually understood in terms of phonons 
(density fluctuations), and within the context of bosonization we give an alternative representation 
in terms of fractional states. This allows to make contact with Bethe Ansatz which predicts similar 
fractional states. As an example we study the spinon operator in the absence of spin rotational 
invariance and derive it from first principles: we find that it is not a semion in general; a trial 
Jastrow wavefunction is also given for that spinon state. Our construction of the new spectroscopy 
based on fractional states leads to several new physical insights: in the low-energy limit, we find 
that the Sz = continuum of gapless spin chains is due to pairs of fractional quasiparticle-quasihole 
states which are the ID counterpart of the Laughlin FQHE quasiparticles. The holon operator 
for the Luttinger liquid with spin is also derived. In the presence of a magnetic field, spin-charge 
separation is not realized any longer in a LL: the holon and the spinon are then replaced by new 
fractional states which we are able to describe. 



PACS Numbers: 71.10 Pm , 71.27-Fa 
I. INTRODUCTION. 

A. Motivations of this work. 

One of the most striking property of some strongly 
correlated systems is fractionalization, that is the exis- 
tence of elementary excitations carrying only part of the 
quantum numbers of the constituent particles of the sys- 
tem. The most famous example is probably the charge 
one-third Laughlin quasiparticle, which is the elemen- 
tary excitation of the fractional quantum Hall fluid at 
filling 1/ = 1/3. El Its existence was recentlw confirmed 
in a beautiful set of shot noise experimentsfl. The ear- 
liest example of fractionalization in condensed matter 
physics is however found in one dimension: the exact 
solution of the Hubbard modelEl by Bethe AnsatafI re- 
vealed that the charge and spin of the electron split into 
two excitations with independent dynamics, known as 
the holon and the spinon. Faddeev and Takhtajan later 
showed that the same spinon is also the elementary ex- 
citation of the ID Heisenberg model: the magnon (the 
usual Goldstone boson) is replaced by two spi|nQns gen- 
erating a continuum for AS — 1 excitationsQu. This 
property of the Hubbard model is known as spin-charge 
separation and is generic of so-called Luttinger liquids 
(LL): LL constitute a universality class for gapless one 
dimensional models such as the Heisenberg chain, the 
Hubbard and t ^ J modelsB. Luttinger liquids are non- 
Fermi liquids: Landau quasiparticlegj are not elemen- 
tary excitations of the LL and as a consequence the elec- 
tron Green's function shows no quasiparticle pole (this 
property is true both for the LL with spin and for the 
spinless LL). Haldane, who coined the name of LL, con- 
jectured that ID gapless models would have the same 
low-energy physics as that of the Tomonaga-Lutttnger 
model. For energies smaller than the bandwidthlljil3, 



the latter naO|del is a fixed point of the renormalization 
group (RG)ll3. In ID, bosonization allows to transform 
the Tomonaga-Luttinger model intOj^ gaussian acoustic 
hamiltonian describing free phononsHJ; the considerable 
success and popularity of bosonization stems from the 
fact that all the computations are straightforward be- 
cause the effective hamiltonian is that of a free bosonic 
field. Another perspective on the LL is provided by con- 
formal field theory (CFT) which describes two dimen- 
sional (or 1 + 1) critical theories with conformal invari- 
ance; this has allowed to identify the Luttinger liquid uni- 
versality class as the set of c = 1 CFTsEj, i.e. the set of 
all models which flas¥ under RG towards the gaussian free 
boson hamiltonianllS. CFT has allowed to formalize the 
finitp-aiTjp analygig of Luttingcr liquids first introduced by 
HaldaneEZnijij. In terms of the gaussian hamiltonian, the 
LL theory can be described as a phenomenological theory 
characterized by the following parameters: u which is a 
velocity for collective modes and K which is proportional 
to the compressibility of the system. 

Yet, although the LL description is supposedly quite 
well established through the formalisms of bosonization 
or CFT, and despite the fact that exact solutions (Bethe 
Ansatz) show the existence of fractional states in the 
spectrum of several Luttinger liquids, there exists no sys- 
tematic study of fractional excitations in the LL to the 
best of the authors' knowledge. What's more, in the 
framework of the bosonization formalism, it is sometimes 
stated that the only physically relevant excitations of a 
LL are phonons, since the effective hamiltonian is just 
that of acoustic phonons. As we show below this state- 
ment is incorrect. Conformal field theory is an alterna- 
tive to bosonization which does stress the spectroscopic 
aspects: yet, application to the study of fractional ex- 
citations in a LL has been limited to the spinon in the 
case of SU{2) symmetry, which is the situation relevant 
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for the Heisenberg chainlli. Fractional excitations must 
exist in a Luttinger liquid if Bethe Ansatz is correct but 
as far as the authors are aware, the characterization of 
these very unconventional fractional states, through ei- 
ther bosonization or CFT, is mostly terra incognita as 
the following list of issues may show: 

1) AS = 1 excitations for the Heisenberg chain form a 
continuum of pairs of spinous. When an Ising anisotropy 
is introduced, in the massless regime (with obvious no- 
tations: \Jz\ < Jx (= Jy) ) the continuum still exists 
and evolves smoothly as a function of the anisotropy 
A = JzJ-jJi. The continuum is again ascribed to pairs of 
spinonaHj. It is intuitively clear that these spinous should 
be in some sense deformations of the SU{2) spinon. We 
will derive in this paper creation operators for these non- 
SU{2) spinous using the bosonization method. 

In the case of SU{2) symmetry, trial wavefunctions 
for these spinous can be found by making use of the 
exact solvability of the Haldane-Shastry (HS) chainEil. 
The HS model shares the properties of the Heisenberg 
chain: it is a gapless SU{2) symmetric spin chain with 
a continuum of spinop-|excitations. Its ground state and 
spinon wavefunction£3 are remarkably similar to those 
one would write for bosonic Laughlin states at filling 
1/ = 1/2: 



*gs(a^i, ..,a;7v) = J|(^j - 2;^)^, (1.1) 

^spinon{zo) = ^(^i ~ ^o) ^(^i ~ ^j)^> (1-2) 

(1.3) 



i<j 



2tt 



Zi = expi—Xi, 

where Xi i^the coordinate a spin down, and a;o that of 
the spinonnJ. We will exhibit a similar wavefunctions for 
the spinon in the absence of SU{2) symmetry. SU{2) 
spinous are semionsE3 (anyons with a statistics interme- 
diate between that of fermions and bosons): we will show 
that the statistics is affected when an anisotropy is intro- 
duced. 

A continuTBa is also found by Bethe Ansatz for ASz — 
transitionsEa. Low-lying excitations are described in 
that approach as two-strings states in the string for- 
malism customary to Bethe Ansatz. This description 
is similar to that given for the A5' = 1 continuum of 
the isotropic chain. In the latter case it is quite clear 
that a spin one-half should be ascribed to each of the 
(pseudo) "hole" states in the string, which leads to the 
spinon interpretation since each state should contribute 
symmetrically to the spin-flipQ. For ASz = transitions, 
the total z spin components of the excitations add up 
to and the continuum results from the excitation of 
particle-hole pairs. For the isotropic chain, this contin- 
uum is generated by spin 1/2 spinon-antispinon pairs. 
By contrast, for the XY chain the ASz = continuum is 
due to particle-hole pairs of magnon-like spin Sz ~ ^ ex- 
citations. The case of the isotropic Heisenberg chain for 
which spin 1/2 spinous are involved both in the ASz = 



and the ASz = 1 continuum, is therefore incidental. The 
important lesson to be learned is that in the presence of 
an Ising anisotropy, the A^z = 1 and ASz — continua 
may involve differen|t-|fractional spin states: in the first 
case we have spinou^3, but in the second case the spinon 
identification is not always correct. What happens in the 
case of an arbitrary anisotropy will be dealt with in this 
paper. 

2) The holon appearing in the exact solution of the 
Hubbard model is a spinless charge one excitatiouQ. 
The issues raised for the spinon (operator, wavefunction, 
statistics) extend naturally to the holon. 

3) Spin-charge separation is an asymptotic property of 
the Hubbard model valid in the low-energy limit. When 
a magnetic field is applied Frahm and Korepin found that 
spin-charge decoupling was not realized even in the low- 
energy limited. In this paper we derive the new excita- 
tions replacing the holon and the spinon. 

4) Xhj the context of the Calogero-Sutherland (CS) 
modelEj the existence of fractional excitations similar to 
Laughlin quasiparticles was suggestedtZl. In a variant of 
the standard LL known as the chiral LL used to describe 
edges of a FQHE sample, Laughlin quasiparticles do ap- 
pear but the existence of such states follows from that of 
the same excitations in the bulkcS. In the CS model the 
proposal was triggered by the similarity of the ground 
state with that of the 2D Laughlin wavefunctions and by 
special selection rules. The ground state iscj: 



'f(xi,..,XN)=Y[\ 



Z{ Zj I , 



1<J 



2n 

Zi = expi—Xi 



(1.4) 
(1.5) 



where A is a coupling constant for the l/r^ intecaction 
potential of the CS model. The CS model is a LlJl3 and 
the LL parameter is just K = 1/A. A pseudo-particle 
formalism similar to that of Bethe Ansatz can be intro- 
duced and for the restricted case of rational couplings 
X ~ p/q special selection rules are found for the dynam- 
ical structure factor: p pseudo-holes (a pseudo- hole is 
a hole in the Fermi sea of pseudo-momenta) must Jae 
accompanied in any excitation by q pseudo-particlestJ. 
For a charge (—1) pseudo-particle, this means one has 
a charge 1/A = K for the pseudo-hole. In the interpre- 
tation of those selection rules it is proposed to view the 
CS model as a gas of non-interacting pseudo-particles 
with anyonic statistics ttA and one rewrites the ground 
state as an anyonic wavefunction Y[i<j i^i ~ ^j)^- The 
pseudo-holes are particle-hole conjugate of these anyons: 
the main modification with the non-interacting case be- 
ing the new selection rule; a wavefunction for these 
pseudo-holes consistent with those interpretations is then 
{zi — z) Y\i<j i^i ~ ^j)^ which has the correct charge 
and statistics. The pseudo-hole is therefore identified as 
a Laughlin quasiparticle. It exists for rational couplings 
and carries the rational charge 1/A. 
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There are however several hmitations to those views. 
Firstly, these considerations are only valid for rational 
couplings (the pseudo-particle selection rules can not be 
extended to irrational A): the physics of the CS model 
is by contrast completely continuous with the coupling 
and does not discriminate between rational and irrational 
couplings. The impossibility to describe rational and ir- 
rational couplings on the same footing means the repre- 
sentation is not adequate. Secondly a disymmetry be- 
tween particle and hole excitations is introduced. 

In a parallel strand of ideas, Laughlin quasiparticles 
were also proposed in studies of transport in a LL. The 
basis of the argument is that for a LL with an impurity 
potential a charge K and not a charge unity is backscat- 
tered at the impurity location (where M is the LL pa- 
rameter, i.e. the conductance of the LL)tJ. The impurity 
potential can be rewritten as a hopping potential for a 
charge K state whose exchange statistics is irK as seen 
from the commutation relations (i.e. anyonic). But for 
K = l/{2n + 1) these states are just those given by Wen 
for the Laughlin quasiparticles of his chiral LL: this sug- 
gests to identify these states as Laughlin quasiparticles. 
The main difficulty in that argument is that it relies on 
the introduction of an impurity potential in the LL: this 
obscures the question of the existence or not of a Laughlin 
quasiparticle in the pure non-chiral Luttinger liquid. In 
summary what is missing is a proof that states similar to 
Laughlin quasiparticles might be exact eigenstates of the 
LL boson hamiltonian (i.e. of the RG fixed point in the 
low-energy limit). The existence of Laughlin quasipar- 
ticles for the non-chiral LL for arbitrary couplings must 
then be considered at this point as an unproved conjec- 
ture. 

The long list of issues we have brought up in points 
(1 — 4) above should convince the reader that a thorough 
discussion of fractional excitations for Luttinger liquids 
within the formalisms of bosonization or CFT remains 
to be done. This is what motivated us to re-examine in 
that paper the spectrum of Luttinger liquids. We want 
to stress that although the previous examples concern in- 
tegrable models, the detailed physics of such integrable 
models is not really our main interest: what matters for 
us is the universal low-energy content of these theories 
and of course we will be unable to tell anything through 
bosonization on the high-energy physics. Although it 
seems to be taken for granted that the excitations of a 
LL are the holon and the spinon on account of Bethe 
Ansatz studies of the Hubbard model, we are not aware 
of any existence proof of such fractional excitations when- 
ever the model is non-integrablc: this is so because any 
proof must resort to the universality hypothesis, that is 
to the LL and bosonization frameworks. The theoreti- 
cal formalism we wish to introduce aims at bridging that 
gap by focusing on the universal structure of fractional 
excitations of Luttinger liquids through the bosonization 
method. 

The structure of the paper will be as follows: section 
m is an introduction to the topics considered in the pa- 



per. In subsection I B we give a short review of the LL 
physics in order to set the notations used throughout the 
paper; the issues discussed in t he pr esent subsection lA 



will be amplified in subsection IB 2 in which we present 
the standard view on excitations of the LL. In the next 
section || we will show that an alternative eigenstate ba- 
sis can be built: that quasiparticle basis allows a natural 
discussion of fractional states. In section [II, we will gen- 
eralize our analysis to the LL with spin. When a mag- 
netic field is added on to the-Hubbard model, spin-charge 
separation no longer occursE3. The standard spin-charge 
separated Luttinger liquid theory is not applicable any 
more. We will introduce in subsection [II B a general 
framework related to the K topological matrix of Wen's 
chiral Luttinger liquidfj, which yields simple criteria of 
spin-charge separation in terms of a Z2 symmetry: we 



will be able then in subsection III C to derive the frac- 



tional excitations which replace the holon and the spinon. 
The general LL theory we have introduced will then be 
applied to th e Hubbard model in a magnetic field in sub- 
section HI D in which we explain the relation between 



our approach and the formalism of the dressed charge 
matrix due to Frahm and Korepin. Let us mention that 
section g of this paper expands on a short yefsion which 
contained results in the case of a spinless LLE3, whereas 
part ^ presents totally new material. 



B. The Luttinger liquid. 

1. Notations. 

This section will define the notations employed 
throughout the paper. We exclusively deal with Lut- 
tinger liquids and therefore when considering some spe- 
cific models such as the Heisenberg spin chain or the Hub- 
bard model we implicitly assume that we are working in 
the LL part of their phase diagrams. The whole physics 
of the LL is embodied in the following hamiltonian: 



Hb = ^ dx K-\V<^{x)f + K{yQ{x)f (1, 
2 Jo 



6) 



supplemented by the so-called bosonization formulas. We 
work on a ring of length L. u and K are the LL param- 
eters. $ can be interpreted as a displacement field for 
phonons, while is a superfluid phase; indeed the parti- 
cle and current densities are defined as: 



p{x) -po = ^V$(a; 

j{x) = -^WQ{x). 
a/tt 



(L7) 
(1-8) 



Renormalized current: Actually j{x) is a bare current 
density which corresponds to the correct one only in the 
non-interacting case K = 1: the continuity equation 
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shows that the correct current density is renormalized 
and is 



jR{x) = uKj{x) 



(1.9) 



(the Fermi velocity has been set to unity) . We will discuss 
in section H the meaning of such a renormalization. 

The particle operators for bosons and right and left 
moving fermions respectively are given as: 

^'5(2:) = : expiV7re(x) :, (1.10) 
^ f,r{x) = : expiv^ (6(a;) — <i>(a;)) : expikpx, (1.11) 
^f,l{x) = ■■ expi^/n {Q{x) + <i>(x)) : exp-ikpx (1.12) 

(in the following we will assume that these operators are 
normal ordered), kp = ttNq/L is the Fermi momentum 
where iVo is the number of particles which is fixed by 
the chemical potential. 8 and 11 = V<i> are canonical 
conjugate boson fields: 



[e(x),u{y)]^i^{x~y)- 



(1.13) 



The zero modes of the charge and current density are 
respectively: 

N^No + Q= [ p{x)dx^No- [ -^V^dx, (1.14) 
Jo Ja V^^ 

J = [ j{x)dx^ [ -^\/e{x)dx. (1.15) 
Jo Jo 



Q has integral eigenvalues as befits a charge operator; 
in the bosonization mapping, the charge quantization is 
taken into account by the topological quantization of the 
phase field $. Similarly, since j{x)dx is a closed line 
integral (around the LL), it is a quantized number: this is 
just the topological quantization of the superfluid phase; 
the normalization of the fields have been chosen so that 
J is an integer. For fermions, Q — N+ + iV_ and J = 
where and iV_ are respectively the (integral) 
number of (bare) electrons added to the ground state at 
the right and left Fermi points. Thciconstruction we have 
reviewed above is due to HaldaneLj. 

Integrating the Fourier expansions of the charge and 
current density gives: 

e(x) = Go + ^Ja; + ^ Ve„expi^a;, (1.16) 
V L — ' '■ 



n^O 



^(x) = $0 - 4=y"*«expi^a;. (1.17) 

^ n^O ^ 

Note that these fields are not periodic: this allows for 
the above mentioned topological excitations. We demand 
that the boson or fermion operators are physical objects 
and be periodic on the ring: "^Bfp{x) = b/f{x + L); 
this then implies the following selection rules on the 
eigenvalues Q and J of the zero modes Q and J: 



Bosons : J even integer, (1-18) 
Fermions : Q — J even integer. (1-19) 

Both Q and J are integers. The zero modes are some- 
times extracted from the definition of the fermion opera- 
tor which defines the U± operators, first built by Heide»|- 
reich and Haldane for the Tomonaga-Luttinger modeltJfl: 



U± = exp iy/n (60 ± ^'o) • 



(1.20) 



It will be useful to consider the commutation properties 
of the following operators: 

Va,f3ix) =: exp -iy/^ (ae(x) - /3$(x)) : . (1.21) 

Using Campbell-Haussdorf formula, one finds: 

Va.^l3{x)Va,i3{y) 



(1.22) 



where sgn{x) is the sign function, which shows in partic- 
ular that ^'_F(a;) is a fermionic operator. We define the 
exchange statistics of an operator per: 

0(x)0(y) 

^ O{y)O{x)exp-i0 sgn{y - x). (1.23) 
For instance 9 — t: for fermions. 



2. Excitations. 

Until Ihe work of Heidenreich0 and subsequently of 
HaldaneO, the only excitations considered in the. gaussian 
model were the bosonic phonon (or plasmonEj) modes. 
But the hamiltonian contains a second part correspond- 
ing to the energies of states with non-zero charge or cur- 
rent with respect to the ground state. In reciprocal space, 
the gaussian hamiltonian becomes: 



^5 = ?E^^"'n,n_ 



g#0 



9l+KJ- 
2L\ K 



(1.24) 



We have split the hamiltonian into the phonon part and 
the non-bosonic zero mode part. The first term can in- 
deed be rewritten as: 



with the phonon operators: 

b„ = 



2 xkr^" 



(1.26) 
(1.27) 
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The second term in the hamiltonian is standard in confor- 
mal field theory; it corresponds to finite size corrections 
to the energy when one adds particles or creates persis- 
tent currents in the Luttinger liquid. The correspond- 
ing states are built by means of Haldane's-I/± operators 
which act as ladder operators in Fock spaceo. This {Q, J) 
part of the hamiltonian is often called in CFT a zero 
mode part. The corresponding excitations may however 
carry momentum. A non zero J excitation creates indeed 
a persistent current with momentum Jkp. These states 
are therefore non-dispersive since their momentum may 
only assume the discrete values Jkp. 

The spectrum of the hamiltonian results from a convo- 
lution of plasmon excitations and of these {Q, J) excita- 
tions as is apparent in figures (0J^): two linear plasmon 
branches rise from each local minimum of the energy ob- 
tained for the zero-mode states (Q, J). It is important 
to note that there are selection rules on the allowed val- 
ues of (Q, J), which refer back to the quantum statistics 
of the particles: as reviewed in the previous section, the 
gaussian model can be considered either for bosons or 
fermions, which results in different bosonization formu- 
las. For bosons, J is constrained to be an even integer 
while for fermions, Q and J must have the same par- 
ity. This then leads to two different spectra as can be 
seen from figures (]l|) and (||): for instance, for bosons 
the state (Q = 1, J = 0) is available while it is forbidden 
for fermions; conversely (Q = 1, J = 1) is available to 
fermions but not to bosons. Thus we have two different 
theories: the same hamiltonian leads to different proper- 
ties depending on whether we considjCr a Fock space of 
bosons or a Fock space of fermion£j. We will call the 
LL with bosonic (resp. fermionic) selection rules: the 
bosonic (resp. fermionic) LL. For the bosonic LL, as de- 
picted in figure (|l|) the spectrum in arbitrary charge sec- 
tors has the same form but for a shift in energies: in the 
charge sector Q one must add the constant ttuQ"^ / {2L) to 
the energy. The same energies are found for the fermionic 
LL in charge sectors for which Q is an even integer, but 
if Q is an odd integer there is a new spectrum with local 
minima at momenta zLkp and not k — (figure |^). 

In the rest of the paper we refer to this parametrization 
of the spectrum in terms of phonons and zero modes as 
the zero mode basis; this is to be distinguished from the 
quasiparticle basis which we will build later. A property 
which will prove crucial for the rest of the discussion is 
the fact that in the free-fermion case a quasiparticle basis 
exists as an alternative to the zero mode basis: instead of 
the zero modes basis, it is indeed possible to parametrize 
the spectrum in terms of the usual Landau quasiparti- 
cles. Below we show that a similar quasiparticle basis 
can be built in the interacting case. While fractional 
quasiparticles do occur in exactly solvable models (the 
holon, the spinon), scant contact had been made with 
the bosonization approach as mentioned earlier. In the 
low-energy limit, using the bosonization formalism, we 
will directly recover the fractional excitations predicted 
in Bethe Ansatz, with the advantage that the simplifica- 



tions brought by the low-energy limit will allow a com- 
plete characterization, giving for instance easy access to 
statistical phases. 



II. FRACTIONAL EXCITATIONS OF THE 
SPINLESS LUTTINGER LIQUID. 

This section is divided as follows: first, we discuss 
the property of chiral separation which is central to the 
physics of fractionalization; then, we exhibit fractional 
quasiparticles for the bosonic LL before turning to the 
fermionic LL for which we will find a different set of ele- 
mentary excitations. 



A. Chiral separation. 



1. Chiral vertex operators and fractionalization. 

The gaussian model is endowed with a very basic prop- 
erty which is that of chiral separation, i.e. we can split it 
into two commuting parts corresponding to right or left 
propagation of the fields. This is a property which is sys- 
tematically used by CFT in the analysis of conformally 
invariant systems. Indeed: 

Hb = ^ dx K-^U{x)^ + i^(Ve(.T))2, (2.1) 
2 Jo 



Q{x,t) = 



(2.2) 



More precisely we introduce the following chiral fields: 

$(x) 



e±(.T) = e(x)T 



K 



which are related to the phonon operators by: 



g > : bq^ 



K\q\ 



e 



+ ,9> 



q<0: 6, = y:^e_,,. 

In terms of these fields the hamiltonian becomes: 
Hb = + iJ_ , 
H±^^ t dx :(9,e±(x))': 

^ \ \ u+h , I Q±Kj\ 



±g>0 



(2.3) 

(2.4) 
(2.5) 

(2.6) 
(2.7) 

(2.8) 



only contains right-moving phonons and similarly 
for iJ_ with left-moving phonons. It is clear also that 
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[H+, H-] = 0. Let us show now that these fields 8± are 
chiral; they obey the equal-time commutation relations: 



e±(x),Ty5,e±(j/) 



i6{x-y), 



(2.9) 



which implies that the momentum canonically conjugate 
to <d± is Ile± = T^dx'd±- The equations of motions for 
these fields are: 



(2.10) 



Thus: &±{x,t) = Q±{x ^ ut) which means that 
we have chiral fields indeed. The superfluid phase 
has therefore been parametrized as: Q{x,t) — 
\ [6+ (x-ut) + Q^{x + ut)] . 

One may define chiral density operators as well as the 
corresponding chiral charges as: 



P±{x) 



2^ 
Q±KJ 



d^^±{x) 



Sp{x) ± Kj{x) 



(2.11) 
(2.12) 



Those chiral densities obey the anomalous (Kac-Moody) 
commutation relations: 

[p±{x),p±{y)]^T-§^dJ{x-y). (2.13) 

Let us now consider the injection of Q particles with a 
momentum q and current J . In that case, the plasmon 
total momentum is equal to q — J{kp + ^). In the 
bosonization formalism, the operator creating this state 
is: 

VQ,j{q) = ^ r rfa:e^(«-^^^)" : exp-iV^{QQ - J$) : . 
VL Jo 

(2.14) 

This can also be rewritten as: 

VQ„j{q)^^ [ da;e^('-'^^^)"exp-iV^g+e+(a;) 
vL Jo 

xexp-i^/^Q-e-{x). (2.15) 
As a function of time: 

VQ,j{q,t) = ^ [ dxe*(«-'^'=^)"exp-^^/^g+e+(.T-M^) 
vL Jo 

X exp —iy/TTQ^Q^{x + ut). (2-16) 

There is therefore a splitting into two counter- 
propagating states. For non-interacting electrons the chi- 
ral charges Q± are integers since K = 1 and the opera- 
tors exp —i^/^^Q±<^±{x^ut) are just those of Q± Landau 
quasiparticless. But in the general case this is not true 
anymore: we will therefore have states carrying fractional 
charges. 



We now define the chiral vertex operators which ap- 
peared in the previous expression as: 



^0+(^) ^ exp-jV7rQ±e±(x), 



(2.17) 



where the upperscript ± refers to the direction of propa- 
gation. They obey the following commutation rules: 



p{x) 



= Q±S{x-y)V^(x), 



Q,V^Jx) ^Q±V^Jx), 



Q± 

K 



(2.18) 
(2.19) 
(2.20) 



which shows they carry charges Q± = Q±^j ^y^i^Yy are 
non-integral in general. The above operator identity 
means that the charge is 'sharp': by 'sharp' we mean that 
the charge found is not a quantum average ( < Q > is not 
necessarily quantized of course). This is a point we want 
to stress because this means that these quantum states 
are genuinely fractional. This shows then that if one in- 
jects Q particles with current J in a LL, one should ob- 
serve a charge (5+ = '^^^■^ state propagating to the right 
at velocity u and a charge Q- = ^^^'^ going to the left 
with velocity ~u. For instance, let us inject an electron 
exactly at the right Fermi point: this is a (Q = 1, J = 1) 
excitation (with no plasmon excited); there would then 
be fractionalization into a charge state going to the 
right and a charge going to the left. 

The most important property of these fractional states 
is that they are exact eigenstates of the gaussian hamil- 
tonian. The proof requires a proper definition of their 
Fourier transform because they are anyons, as will be 
shown shortly: from equation ( |l.22| ) it is clear indeed 
that the commutation relations are anyonic with an any- 
onic phase 



= ±11- 



K 



(2.21) 



Due to its anyonic character Vq^{x) does not obey pe- 
riodic boundary conditions; if we use t he expressions of 
the fields $ and 9 ( equations |1.17| , |1.16 ), we immediately 
find that: 



V^^ {x + L) = exp±z27r^ V^^ (x). 
The Fourier transform is then defined as: 



with a pseudo- momentum qn quantized as: 
27r 27r (3± 

^'^-^TIT' 



(2.22) 



xV^Jx), 



(2.23) 

(2.24) 
(2.25) 
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(where we have defined a phonon part Tj^ of the momen- 
tum). 

The operators Vg^^(g„) are such that: 
1) VqIj^ (g) 1 5*0 > is an exact eigenstate of the chiral 
hamiltonian H± with energy: 



EiQ±,qn) 



(2.26) 



where \'^o > is the interacting ground state (see ap- 
pendix). It has a linear dispersion. 

2) The states created by the VJ^^((7„) to which one 
adds the phonon excitations form a complete set. This 
is obvious because the states Vq^j{x) span the full Fock 
spaceEa. 



2. The LL spectrum in terms of fractional quasiparticles. 



Let us consider figures (|l|,y) which show the spectrum 
of the LL hamiltonian in various charge sectors and ask 
the following question: what happens when one adds Q 
particles to the system (i.e. in the charge sector Q)7 
In the standard view of the LL spectrum based on the 
phonon and zero modes basis, the dynamics of the charge 
added to the LL is unclear because it is concealed in 
the zero modes. The parametrization of the spectrum in 
terms of the zero modes and the phonons does not allow 
to find what happens once the charge Q is added to the 
system because that choice of basis involves the use of 
non-dynamical states (Haldane's Upm operators, which 
describe the zero modes). By contrast the quasiparticle 
basis only involves states which have a dynamics (the 
phonons and the fractional states) and we are therefore 
able to tell what happens to the charge, how much of it 
will move to the right, and so forth: if we consider the 
two branches starting from k = Jkp in the charge sector 
Q, the right branch corresponds to a right-moving frac- 
tional excitation with linear dispersion and with charge 
Q+ — ^^Y^j while the left branch is due to a left-moving 
fractional state with charge Q- ~ ^^^^ ■ The contin- 
uum in between the branches simply results from the cre- 
ation of the two fractional excitations with both non-zero 
momentum q+^JT and q-.n (on the right branch, a charge 



Q- = 



}-K.l 



is also created but it has zero momentum 



g__„ ~ 0, and conversely on the left branch). 

The direct way to find out how the charge Q will be- 
have, is to exhibit the quantum states which will describe 
the propagation of the charge. This is what the quasi- 
particle basis does because it directly considers the states 
involved in the dynamics of the charge. Of course, the 
two bases (the quasiparticle basis and the zero mode ba- 
sis) are mathematically equivalent and therefore lead to 
identical physics: therefore the charge dynamics can also 
in principle be determined in the zero mode basis, but in 
the quasiparticle basis, we have the benefit that the spec- 
troscopy immediately tells us the fate of the charge added 



to the system. In sharp contrast, in the zero mode basis, 
the spectroscopy is not useful because the states used in 
that basis are the phonons (which have no charge) and 
the U± operators (which have charge but no dynamics). 
We will give such an argument in the next section: this 
will prove in an independent manner the fractionalization 
of the LL spectrum (in a way which does not depend on 
the explicit construction of the fractional states opera- 
tors). 



3. Selection rules and fractionalization. 

The fractional charges carried by the fractional excita- 
tions considered above are not arbitrary: they must take 
on the values 



Q± = 



'±KJ 



(2.27) 



where both Q and J are integers. We may view these 
constraints on the allowed spectrum of fractional charges 
as selection rules. These selection rules have however a 
clear physical meaning which we discuss now. 

Although these excitations do not carry the electron 
quantum numbers because of the fractionalization of the 
spectrum, nevertheless the elementary constituents of our 
system are electrons (they are the high energy elemen- 
tary particles of our systems): this means that they alone 
define the structure of Fock space, with the implication 
that all physical states must consist of an integrer num- 
ber of electrons. Despite the fact that there are fractional 
states, the previous remark implies that these fractional 
states will be created in appropriate combinations so that 
the total charge is always an integer. This is the expla- 
nation of the previous selection rules we found, which in 
fine enforce the basic constraint that we started out with 
electrons. We may view these selection rules as being 
topological since they are directly related to the struc- 
ture of the Fock space. 

It is easy to show that eq. ( 2.27 ) immediately follows 
from the requirement that all states are electronic. We 
consider two counterpropagating states with arbitrary 
charge and Q-; we make no hypothesis on the val- 
ues of the charges, nor on the nature of the chiral states 
(we do not assume they correspond to Vq^ and Vq ). 
The only assumptions we make are the following: (a) 
the one-dimensionality which means that the eigenstates 
have momenta in one of either two directions and (b) 
that the current density operator is renormalized. We 
then have two constraints on the values that the charges 
(5+ and Q- may assume: since our Fock space is that 
of electrons, all the states contain an integer number of 
electrons i.e. Q+ + Q_ = Q is an integer. The second 
constraint stems from the renormalization of the current 
density operator: 



jR{x) = uKj{x) = uK 



1 



(2.28) 
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where j{x) is the current density in the non- interacting 
case. This expression can be derived from the continuity 
equationla. Going around a ring of length L in the LL we 
get a (persistent) current which must be quantized: 



Jr 



dxjji{x) = uKJ 



(2.29) 



where J is an intege but the current carried by the 
states with charges (5+ and Q_ is Jr = u{Q+ — Q-). 
Therefore: 



while: 



{Q+ — Q-) — KJ, J integer 



(Q+ + Q-) = Q, Q integer 



(2.30) 



(2.31) 



Solvin g for these constraints, one recovers the selection 
rules (2.27), i.e the spectrum of fractional charges. We 
observe in passing that this argument does not depen d on 
our formal algebraic derivation of subsection (II A 1) and 
provides an alternative proof of the existence of fractional 
states as well as it yields the allowed charge spectrum. In 
that argument, fractionalization follows from the renor- 
malization of the current in the presence of interactions. 



B. Elementary excitations of the bosonic LL. 

1. Elementary excitations. 

We now establish a series of new results concerning 
the elementary chiral excitations of a non-chiral LL. We 
would like to find a basis of elementary excitations, i.e. 
identify objects from which all the other excitations can 
be built. It will be useful to use a spinor notation to 
represent the fractional states: 















)- 













(2.32) 



(2.32) should be understood as follows: the fractional 
state Vq^ which is an anyon propagating with velocity 

u is created along with the fractional state Vq which 
propagates in the opposite direction with the velocity —u. 
The selection rules are encoded in the second spinor: the 
equation is then read as meaning that addition of Q par- 
ticles with (persistent) current J will result in a splitting 
into the two counterpropagating fractional states Vq"^ 

and Vq_ . 

We must carefully distinguish between Bose and Fermi 
statistics because of the constraints on Q and J. Let us 
consider bosons first: since J is even we can rewrite it as 
J = 2n where n is now an arbitrary integer. But then 
for bosons this implies that the spinor can be written in 
terms of two other independent spinors: 



Q+ 
Q- 



■ Q+KJ 

Q-KJ 
2 

.2 



K 



(2.33) 



This implies that in real space the fractional charge 
excitation is: 



(2.34) 



where Q and n are now independent integers of arbitrary 
sign: {Q,n) £ Z'^. In reciprocal space, one has a convo- 
lution for the exact fractional eigenstate: 



VqA<1) = 



Q 



n [v,%m\ n i^iKipj)] 



Q n 
i=l j=l 



(2.35) 



(the momenta in that expression are the phonon parts 
of the momentum of the operator: for Vq^ {in), In = 

^±^^andg;r=^)i. 

The above equation demonstrates clearly that the ex- 



citation Vq (q) can be built from Q charge 1/2 states 



V^2 ^-nd n charge zLK states V^j^. The whole spectrum 
of fractional excitations is thus built by repeated creation 
of Euid V^j^ which means that they are the elemen- 
tary excitations we were seeking. These two elementary 
excitations will be identified in the following as respec- 
tively the spinon (for spin systems) and a (ID) Laughlin 
quasiparticle. 



2. Wavef unctions of the fractional excitations. 

To be complete, we compute the wavefunctions of all 
the chiral excitations. We will first need the ground state 
wavefunction which is simply a Jastrow wavefunction: 
this is of course expected since the gaussian hamiltonian 
is the ID version of the acoustic haipiJtonian of Chester 
and Reatto's Jastrow theory of He4c3 and is also iden- 
tical tOiBohm-Pines RPA plasmon hamiltonian adapted 
to 1DE3. Since the gaussian hamiltonian is a sum of os- 
cillators, the ground state is a gaussian function of the 
densities: 

*0,s({Pg}) 

= exp(- J2 ^^iP-jP-i) (2-36) 

= exp — — / / dxdx' p{x) In sm—{x — x') p{x'). (2.37) 
J J E 

This expression is valid for the bosonic LL; for the 
fermionic LL, antisymmetry is recovered by observing 
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that the fermionic LL simply derives from the bosonic LL 
through a singular gauge transformation [©n the bosonic 
LL (the Jordan- Wigner transformation)L3: this is ex- 
actly as in the composite boson Chern-Simons theory 
for which the hamiltonian is plasmon-like at the one- 
loop level (RPA) and whose ground state is of course 
symmetric (the modulus of Laughlin wavefunction) ; in 
that theory, the Laughlin state is then found-after undo- 
ing the Chern-Simons gauge transformationEJ. Similarly 
undoing the Jordan- Wigner transformation amounts to 
multiplying the bosonic ground state by the phase fac- 
tor U^<j sgn{x, - Xj) = n,<j (it^) (tliat phase fac- 
tor is found by applying the Jordan- Wigner operator on 
the ground state). The wavefunction is the ID analog 
of the 2D Laughlin state of FQHE if we rewrite the 
previous expression in terms of the particles' positions: 
p{x) = 8(x — Xi), and by introducing the circular co- 
ordinates z = expi^x: 



iPo,b{{xi}) = Yl 



,1/K 



(2.38) 



The wavefuactions of the excited states can now 
be computedlid. Let us consider first the operator 
e'xp —iy/TraQ{xo); since 6 is the canonical conjugate of 
the field U^d^<^> ^ -\^Sp, 



Six) 



1 



iy/n Sp{x) ' 



(2.39) 



and therefore: 



exp —^^/^^aQ{xo)\'^o,B > 



expa 



6p{xq) 



xexp^ p»ln|sin^(.-x')p>') 



K 



exp — / da; p(a;) In sin —{x — xq) 



^0, 



,a/K 



Y[\z^-Zj\ 

i<j 



1/K 



(2.40) 



where C is an unessential constant. 
Similarly: 



.Q± 



eyip±iy/n——^{x) expq=i— — fc^cc 
K K 

expTiTT^ J p{y)dy 
expT^TT-^ j p{y)6{x-y)dy 



n 



{Xi - x) 



\xi - x\ 



(2.41) 



n 



(2.42) 



where we use the definitions of p and z introduced above 
and where 



0{x) 



1 



■In 



(2.43) 



is the Heaviside step function. 

The above operator can thus be seen as a general- 
ized Jordan- Wigner operator, since it multiplies wave- 
functions by a singular phase factor; in this manner we 

recover the phase Y[i<j j^'-z^ l ) the ground state of 
the fermionic LL. Finally we have that: 

^Q+(^)*o,B(a;i, ..,XAr) 



X expikp 



0+ 

2K 



No 



(2.44) 



with a similar expression for Vq (the bar over z denotes 
complex conjugation). It is noteworthy that these wave- 
functions are obtained by multiplying a Jastrow ground 
state with a Laughlin- like prefactor Yii \zi — z^^^^ 
which generalizes the Laughlin quasihole factor OiC-^i ~ 
z). We can now write down the wavefunctions of the two 
elementary excitations: 



^i/2(2^)*o,B(a;i, ..,xjv) 



X exp —I 



ni 



i<j 



2K V N> 



(2.45) 



and. 



V^{x)^a,B{xi, ..,xn) 

c]i{z,^z) Hi 



X exp —ikp 



Zl 
i<j 

No 



,1/K 



(2.46) 



We see that V^{x) is nothing but the ID counterpart 
of the 2D Laughlin quasi-hole wavefunction, provided 
we make the following correspondence between ID and 
2D wavefunctions: K <;=4> z — expi27rx/L <^=^ 
z — x + iy (up to a galilean boost absorbing the factor 



exp —ikp 



No 



X \): in view of the formal analogy 



we will call that state a ID Laughlin state. 
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3. The spinon. 

We found an elementary excitation Vy^{x) for the 
bosonic LL carrying a charge 1/2. When we consider 
spins, which are hard-core bosons, this result translates 
into having a state carrying a spin ASz — 1/2 with re- 
spect to the ground state. In spin language, adding a 
particle into the system {Q = 1) corresponds to fli pping 
a spin {ASz = 1). But it follows from eq.(|l33|) that 
this excitation is a composite of two elementary exci- 
tations, each carrying a charge 1/2. Therefore a pair 
of states with spin Sz — 1/2 is created when one flips 
a spin {ASz = 1)- We naturally identify this fractional 
spin excitation as a spinon. The spinon can be generated 
without any Laughlin quasiparticless if the spin current 
is zero (J = 0): this is a process which we term a pure 
spin process, to be distinguished from a pure spin current 
process {Sz = 0) which generates Laughlin quasiparticle- 
quasihole pairs (see below). 

The properties of the spinon specifically depend on the 
LL parameter K. Although the spin is always Sz — 1/2 
the exchange statistics varies continuously with K (i.e. 
when one varies the interaction): 



4:K 



(2.47) 



For instance for K = 1/2 (which corresponds to SU{2) 
symmetric spin interactions) the spinon is a semion. In 
that special case, the spinon wavefunction we obtain co- 
incides exactly with that ptoposed by Haldane for the 
Haldane-Shastry spin chained: 



xexp-i/cj. (^^^H-x^ (2.48) 

In this expression the coordinates are those of the down 
spins. For K — 1/2 the spinon and the spin K Laughlin 
quasiparticles are identical. Although we have discussed 
fractional excitations for spin systems, the previous con- 
siderations apply of course to bosons: the "spinon" is 
then a charge 1/2 excitation. For convenience we will call 
the excitation a spinon even when we consider bosons. 



Jf.. The LL Laughlin quasiparticle. 

The second elementary excitation we found has the 
following wavefunction: 



Laughlin— qp 



{z,) = X{{z,-z,) Hi 



,1/K 



X exp —ikp 



i<j 
Nq 



+ xo], (2.49) 



which leads us to identify it with a Laughlin quasiparti- 
cles. The parallels which can be drawn between the 2D 



Laughlin quasi-hole and the Luttinger liquid Laughlin 
quasiparticles are indeed very strong. For instance as in 
2D one can use the plasma analogy to find the fractional 
charge: 



l"^ Laughlin~qp{zo)\ — 

exp j j dxdx' [p{x) + KS{x — xq)] 



Yl{z^ - Zq)YI\z, - Z. 

i i<j 



,1/K 



In 



sm-{x-x') [p{x') + KS{x' -Xq)]. (2.50) 



The above expression clearly shows that the charge car- 
ried by the excitation is K in agreement with the direct 
algebraic determination (using the operator Vk)- There 
are however several differences between the LL Laughlin 
quasiparticles and its 2D famous counter-part; first, there 
is no analyticity requirement in the ID problem, since we 
do not have to project into the lowest Landau level: we 
have two chiralities and the LL Laughlin quasi-electron 
is simply 



Laughlin-qe{zQ) = I^l^J ~ Zq) ^ J]^ | 



,1/K 



X exp —ikp 



xo . (2.51) 



Second, while topological quantization forces the 2D 
FQHE Laughlin quasiparticles to have a rational charge, 
the charge of the ID LL Laughlin quasiparticles can take 
on arbitrary real positive values, in particular irrational. 
This is a very startling property: irrational spin had 
already been proposed for solitons in coexisting CDW- 
SDW systems by B. Horowitala, but in a sense this is 
perhaps less surprising since in one dimension there is 
no quantization axis for spin which can therefore take a 
continuum of values. We show below that the Laughlin 
quasiparticles also exist for the fermionic LL; furthermore 
we will find that for the fermionic LL there is another ele- 
mentary excitation which may have an irrational charge. 

How are Laughlin quasiparticles created in a LL? They 
are generated whenever J 7^ 0; they are always created 
as quasiparticle-quasihole pairs. In particular in pure 
current processes {Q = 0) no "spinon" is created and we 
have only Laughlin quasiparticle-quasihole (qp-qh) pairs. 
For a persistent current J excitation with Q = it follows 
from the expression Q± = ^^^-^ that J/2 quasiparticle- 
quasihole pairs are generated. 

From the above analysis we now can give a physical in- 
terpretation to the renormalization of the current density 
operator in the presence of interactions: 



= -^dxQ{x), 
uK 

— ^ iR{x) = —^dxQ{x). 



(2.52) 
(2.53) 
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The velocity u has been normahzed to the Fermi velocity 
so that u = 1 in the absence of interactions for fermions 
{K — 1). We have found that current excitations were 
due to Laughhn quasiparticles. The natural explanation 
of the renormalization is therefore that the current is no 
longer carried by Landau quasiparticles but by Laughlin 
quasiparticles with velocity u and charge K. 



5. The hosonic LL spectrum in terms of fractional 
elementary excitations. 

For the bosonic LL we can now add the following pre- 
cisions to the description of the spectrum. 

For Q — G excitations (see figure (|l|)), the continuum 
is due to multiple Laughlin quasiparticle-quasihole pairs: 
the right branch starting at fc = 2A;f corresponds to the 
propagation of a ID Laughlin quasielectron while the left 
branch is due to a Laughlin quasihole; in between the 
two lines, we have a continuum generated by these two 
excitations. More generally at fc = 2nkp where n is an 
arbitrary integer, the two branches create a continuum 
of n Laughlin quasiparticle-quasihole pairs. (Note that 
for n = 0, which is an exceptional case, we have multi- 
phonon processes.) Therefore the spectrum in the zero 
charge sector is not a Landau quasiparticle-quasihole pair 
continuum except at the special value K = I which de- 
scribes indeed in the low energy limit a gas of hard-core 
bosons. 

In the charge sector Q = 1, pairs of charge one-half ex- 
citations are created: they correspond to the "spinons" of 
spin chains; the pairs are superimposed on the previous 
Laughlin quasiparticle continuum: for instance a 2fci? ex- 
citation generates a Laughlin quasiparticle-quasihole pair 
in addition to the "spinon" pair. 

The Laughlin quasiparticle and the spinon are dual 
states for the bosonic LL; by duality we mean electro- 
magnetic duality which exchanges charge and current 
processes. Indeed the "spinon" is associated with charge 
processes while the Laughlin quasiparticles is due to cur- 
rent excitations. The duality operation which maps a 
bosonic LL onto another bosonic LL is:K = 1/{4K') with 
O = 2$' and $ = 0'/2; zero modes then transform as 
J — 2Q' ct Q ~ J'/2. With these relations, the selec- 
tion rule remains bosonic {J' even) while the hamiltonian 
Hb [K, 6, $] = Hb [K' , 6', retains a gaussian form. 
It is clear then that K — \/2 is, & self-dual point while 
Vk and V1/2 create dual quasiparticles. This is not true 
for the fermionic LL. 



6. The XXZ spin chain. 

Let us illustrate these results on the specific example of 
the anisotropic Heisenberg XXZ spin chain. The hamil- 
tonian of the XXZ spin chain with anisotropy A, after 
a bipartite rotation is: 



if [J, A] = |-i (5+5r+i + Srst+,) + ^StStA . 

(2.54) 

As A is varied, one finds three phases: i) for A > 1 
one gets an Ising antifcrromagnet the twofold degenerate 
groundstate of which leads to solitonic excitations with 
spin one-half 1/2 domain walls; ii) for A < 1 one has an 
Ising ferromagnet; iii) for -1 < A < 1 we have the so- 
called XY phase: this is the Luttinger liquid phase we 
are interested in. The isotropic Heisenberg chain with 
SU{2) invariance corresponds to A = 1 . The Luttinger 
liquid parameter was determined exactly by Luther and 
PeschpLon the basis of a comparison with the Baxter 
modeled: 



K{A) 



2 arccos (—A) 



(2.55) 



The spectrum in the sector ASz = 1 is shown in figure 
(^) for the Heisenberg model; its linearization through 
bosonization is also shown in figure. 

Given thaLja spin one-half can be mapped onto a hard- 
core boson:o through the Holstein-Primakov transfor- 
mation, we can transpose the results we found for the 
bosonic LL to the XXZ spin chain. 

If we want to compare the bosonization linearized spec- 
trum to the exact one there are however two provisos: (a) 
we have to shift the bosonization spectrum by a momen- 
tum TT : this is due to the bipartite transformation one 
makes in the bosonization of the XXZ spin chain (in or- 
der to change the sign of the XY term) and (b) , there is 
a Brillouin zone: therefore we have to identify momenta 
modulo 2tt and since the Fermi vector is kp — tt/2, ex- 
citations with J/2 odd (resp. even) correspond to the 
same harmonics k = n + Jkp = (resp. tt). Taking 
(a) and (b) into account , we can use the results of the 
previous section pertaining to the bosonic LL, to recover 
the linearized spectrum of the XXZ chain. 

We first consider the spin sector ASz = 0. In figure 
(^) starting from momentum tt we have two straight lines 
corresponding to left and right moving phonons, bound- 
ing a continuum; due to the folding of the continuum 
spectrum of the bosonic LL, one superimposes on these 
lines the lines due to the creation of any even number of 
Laughlin qp-qh pairs (the qp dispersion being given by 
one line, and that of the qh by the other; if the qp is right- 
handed, its dispersion is that of the right line, etc.). 
Similarly the lines starting from momentum zero or 2tt 
correspond to the creation of an odd number of Laugh- 
lin qp-qh pairs. The continuum is therefore seen to be 
parametrized entirely in terms of the phonons and Laugh- 
lin quasiparticle-quasihole pairs while the zero mode ba- 
sis relies on phonons and zero modes. The ASz = 1 
continuum is described in a similar manner but for the 
substitution of the phonons by a pair of counterpropagat- 
ing spinons. In the special case of SU{2) symmetry, the 
Laughlin quasiparticle and the spinon become identical 
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operators. The previous parametrization reduces then to 
one involving only pairs of spinous because a pair of coun- 
terpropagating spinon plus a pair of counterpropagating 
spinon-antispinon is equivalent to a pair of spinons prop- 
agating in arbitrary directions. One then recovers the 
Bethe Ansatz result. 

In the low- energy limit, we can now answer the various 
questions raised in the introduction about the spectrum of 
the XX Z chain: 

-what is the nature of the ASz = 1 continuum? It 
is indeed a spinon pair continuum; but superimposed 
on them Laughlin quasiparticle-quasihole pairs can ex- 
ist. The spinon changes when the anisotropy is varied: 
it acquires a statistical phase ■n/AK ~ arccos(— A)/2. 
Therefore the spinons at different anisotropy are r^t adi- 
abatically connected: they are orthogonal states^ this 
is consistent with numerical computations of the spectral 
density of the SU{2) spinon, where it is found that the 
SU {2X-spinon has a zero quasiparticles weight for the XY 
chairiLB. 

-what is the nature of the ASz = continuum ? It is a 
Laughlin qp-qh pair continuum with an unquantized spin 
Sz = = ±7r/(2 arccos (— A)); in the SU{2) symmet- 
ric case, they are identical to the spinons. In the XY 
limit, one recovers the standard spin one continuum pre- 
dicted through a Jordan- Wigner transformation (K — 1, 
Sz = But in between these two points, the elemen- 

tary excitation is neither a spinon nor a Jordan- Wigner 
fermion. 



C. The fermionic Luttinger Liquid. 

1. Elementary excitations : the Laughlin quasiparticles and 
the "hybrid state". 

We now turn to fermions; the analysis of the elemen- 
tary excitations will differ from that found for the bosonic 
LL because the allowed {Q, J) states obey different selec- 
tion rules, namely J is not constrained any more to be 
an even integer, but must have the same parity as Q. We 
m ay th erefore write Q — J = 2n. Then for fermions using 
eq. 



2.27 



{Q+,Q-)^Q 



l + K 1-K 



n {K,-K). (2.56) 



The most general excitation once again, is built by ap- 
plying Q times and/or n times V^j^ to the ground 

2 

state ; this means that we have identified a set of elemen- 
tary excitations for the fermionic LL. Here too we find 
Laughlin quasiparticles V^j^, but instead of the spinon 
we get a "hybrid state": this is a consequence of statis- 
tics; as we will show below, that hybrid state is self-dual 
and is intermediate between the Laughlin quasiparticle 
and its dual state. 

The Laughlin quasiparticle is created by current exci- 
tations: for a pure current process (Q = 0, J 7^ 0) one 



indeed generates Laughlin qp-qh pairs as the above equa- 
tion shows. The continuum for zero charge excitations 
{Q — 0) is often depicted as a (Landau) particle-hole con- 
tinuum as in the non- interacting system (K = 1): this 
is incorrect; wc have instead a Laughlin quasiparticle- 
quasihole continuum. The latter does reduce to the 
standard Landau quasiparticle continuum when K — 1. 
For k = Jkp there is a local minimum of the energy 
from which two linear branches rise corresponding to 
J/2 = —n pairs of Laughlin quasiparticles and quasi- 
holes. For the fermionic LL, the Laughlin quasiparticle 
is not the only state which may have an irrational charge: 
this is also possible for the hybrid state. 

The hybrid state is created in mixed charge and current 
processes: this is the main difference with the bosonic 
LL for which the decoupling between charge and current 
processes is complete. As reviewed in the introduction, 
there are even-odd effects in the fermionic spectrum: the 
spectra obtained by adding an even or an odd number of 
particles are qualitatively different for the fermionic LL. 
The Q = 1 continuum is understood as follows: the two 
branches at kp correspond to a pair of hybrid excitations 
carrying a charge and ^-y^, and propagating with 
velocities respectively —it and u . At —kp the correspon- 
dence is reversed. More generally at Jkp (J is an odd 
integer if Q = 1), in addition to the hybrid quasiparticles, 
one also creates J/2 — Q pairs of Laughlin quasiparticle 
and quasihole. When K = 1 the hybrid states reduce 
to Landau quasiparticles. It is interesting to note an 
evidence for these states in the work of Safi and Schulz 
who considered the evolution of a charge 1 wavepacket in- 
jected at kp in a LL: they found that there was a splitting 
with an average charge < Q >= {l-\-K)/2 propagating to 
the right ;uijd an average charge < Q >= {1 — K)/2 going 
to the leftcZl. This is exactly what we predict. Note how- 
ever a crucial difference: the charge they find is a quan- 
tum average while we deal with elementary excitations 
(exact eigenstates); this has an important consequence: 
while it is clear that on average a charge may assume ir- 
rational values, our result goes beyond that observation 
since it proves that there may exist in condensed mat- 
ter systems a genuine good quantum state with sharp 
irrational charge. 

In this section, we have found that for the fermionic LL 
there are two elementary excitations. One is the Laugh- 
lin quasiparticle already found for the bosonic LL. The 
second one is a hybrid state intermediate between the 
spinon and the Laughlin quasiparticle. The excitations 
corresponding to Q = transitions (they are particle-hole 
excitations in the non- interacting case), form a Laughlin 
quasiparticle-quasihole pair continuum when K 1. 



2. Dual basis and the dual quasiparticles. 

The elementary excitations we have derived form a ba- 
sis from which all the LL spectrum is recovered; by no 
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means is this choice of basis unique: other bases of el- 
ementary excitations are generated with matrices asso- 
ciated with a basis change having integer entries whose 
inverses are also integer-valued: this ensures that all ex- 
citations are integral linear combinations of the elemen- 
tary excitations. (The matrices belong to 5i(2, Z).) For 
instance for fermions, another basis of elementary exci- 
tations consists of states and V-^: 

2 

iQ+,Q-) = J (^^, + n (1, 1) . (2.57) 

It is actually a dual basis to the previous one: for 
fermions the electro-magnetic duality which exchanges 
charge and current excitations is expressed by K < — > 
1/K and $ < — > 0. This is a canonical transforma- 
tion; it results in: Hg [K, 6, $] = Hb [1/K, $, 6]. The 
fermionic selection rule Q — J even is obviously preserved: 
the transformation is therefore a duality operation for 
the fermionic LL. Observe that the transformations dif- 
fer from those for the bosonic LL. 

What is the nature of the elementary excitation V-^? 
Under the duality transformation — > V^^ and 
V^j^ — > V^^. Therefore is an excitation dual to 
the Laughlin quasiparticle. It carries a charge unity and 
its wavefunction is: 

v,Hzo)^F = n(-^ - -o)^/^ n - -^i'^" n 

■I 1<J 1<J 

(2.58) 

We stress that although carries a unit charge, it is not 
an electron: the statistical exchange phase is n/K which 
means Vf^ is an anyon (it can be a fermion, in the special 
case K ~ l/(2ri + l)). The difference with the electron 
is quite clear since the electron creation operator is: 

^(x) = ^ expi(2n + l)kpx 

n 

expi (V^e(a;) - V7r(2n-f l)$(a;)) , (2.59) 

while: 

y+(2;) =: expiv^(e(2;) - ^5^) : . (2.60) 

For K 7^ l/(2n + 1) the dual excitation appears to be a 
non-linear soliton of the electron. This excitation is inter- 
esting in many respects, li K = l/(2ri-|- 1) (the Laughlin 
fractions) the excitation is fermionic and the exchange 
statistics of the operator is 7r(2n + 1). The dual quasi- 
particle corresponds then to a sub-dominant harmonic of 
the electron Fourier expansion around k ~ (2n -I- l)kp. 
If one attaches 2n flux tubes to the electron (i.e. multi- 
plies the electron operator by the Jordan- Wigner phase 
exp i^/T:2n<^) the dual state becomes the dominant k — 



kp harmonics: this is exactly the composite fermion con- 
struction and it may then be more fitting to speak of a 
composite fermion (indeed, the statistics of the operator 
is (2n -I- l)7r and not tt). Because of the similar long dis- 
tance behavior of their Green functions. Stone proposed 
to identify such a sub-dominant operator -which he calls 
a hyperfermion-pHfith Wen's electron operator introduced 
in the chiral LLc3. This hyperfermion is identical to the 
dual state for K = l/(2n -f 1). In general the dual state 
and the electron are however orthogonal: this is quite 
clear when one considers the LL with spin. The dual 
state is then generalized to a state with the same quan- 
tum numbers as the electron (carrying a unit charge and 
a spin one- half), but with again anyonic statistics. But 
due to spin-charge separation, that state is not stable 
and decays into a spinless charge one quasiparticle which 
is none other than the holon, and a spin one-half ex- 
citation, which is just the spinon. The dual excitation 
we have found is therefore the analog of the spinon and 
the holon for the spinless LL and has nothing to do (in 
general) with an electron. In the following in accordance 
with the previous remarks we will call these states holons 
(for the spinless LL) or dual states. 

These dual holon states also occur in Haldane's in- 
terpretation of the Calogero-Sutherland model: he pro- 
posed that a natural interpretation of such a model was 
not in terms of elertrops or bosons but as a gas of non- 
interacting anyon£3'tZl. The basis for that interpreta- 
tion is the finding that for rational values of the coupling 
X — p/q {X is related to the LL parameter by the sim- 
ple relation A = l/K), the dynamical structure factor 
obeys simple selection rules: q "charge —1 bare parti- 
cles" (the anyon -it has anyonic statistics vrA — tt/K) 
are created with p " holes" which therefore carry a charge 
1/A. The particles and holes appear as pseudo-particles 
in a pseudo-momenta parametrization of the spectrum. 
There are however difficulties with that interpretation: 
it is not clear how the selection rules are generalized 
when the coupling A is irrational; the physics is indeed 
completely continuous with the coupling while the se- 
lection rules are only valid for rational couplings; be- 
sides, an asymmetry is introduced between particles and 
holes. This means that this parametrization which re- 
lies on pseudoparticles is probably inadequate. In the 
low energy limit, the Calogero-Sutherland model has the 
properties of a Luttinger liquid. It is therefore possi- 
ble to describe its quasiparticle spectrum in terms of the 
fractional excitations we have found in this paper: we 
do find the charge 1 anyon proposed (this is our dual 
state); however our selection rules are quite different. 
First they depend on the statistics (electrons or bosons): 
in contrast the pseudo-particle based selection rules do 
not involve spinons; this runs contrary to results on the 
bosonic LL for which the dual state must actually be 
seen as a composite state made out of two spinons. Sec- 
ond, our selection rules are valid even for irrational cou- 
plings that is imply the existence of quantum states with 
sharp irrational charges. Third our selection rules respect 
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particle-hole duality: there exist both a charge one anyon 
with statistics tt/K (the dual holon) and a similar charge 
— 1 anyon; the same applies to Laughlin quasiparticle for 
which we have both quasiholes and quasielectrons. Our 
selection rules involve the hybrid state and the Laughlin 
quasiparticles, and the holon only appears in the dual 
basis where it is accompanied again by the hybrid quasi- 
particle. 

Finally, a symmetric basis can be associated with the 
hybrid excitation: 



Q+ 
Q- 



l+K 



l-K 

l±K 
2 



(2.61) 



where to = 2~ ^'^'^ ^ ~ 2 ' ^'^^ again in- 

dependent integers. They physically correspond to the 
number of electrons added to the system at the right 
and left Fermi points respectively. That self-dual basis 
reduce to Landau quasiparticles when K = 1. The phys- 
ical processes generated in that symmetric basis are not 
charge or current excitations but addition of electrons 
at the Fermi surface. Note that the arbitrariness in the 
choice of a basis simply reflects the possibility to stress 
various specific physical processes as elementary. But ex- 
periments probe and (5_; for a given set of Q and 
J, Q± assume the same value irrespective of the basis 
choice. 



III. THE LUTTINGER LIQUID WITH SPIN. 

In this section we generalize the construction of frac- 
tional excitations developed in section (|l|) to the full Lut- 
tinger liquid with spin. One of the main properties ex- 
hibited by the effective theory is spin-charge separation, 
the complete decoupling of spin and charge dynamics. 
In the exact solution of the Hubbard model by Bethe 
Ansatz, excitations display such a spin-charge separa- 
tion: one state, the holon is a spinless particle carrying 
the charge of the electron, while the other, the spinon, 
is a neutral spin one-half stateEl. This is an asymptotic 
property only valid in the low-energy limit (the Hub- 
bard model in the large U limit is an exception because 
spin-charge separation is realized at all energy scales). 
However this property is not obtained for all the gapless 
itinerant ID models, in the low energy limit. According 
to the universality hypothesis they should be described 
by the Luttinger liquid framework if the interaction is 
not too long-ranged. One such example is the Hub- 
bard model in a magnetic field which does not display 
spin-charge separation even in the low-energy limit, al- 
though it is a short range gapless model. This model 
was analyzed by Frahm and Korepin in the. framework 
of Bethe Ansatz plus conformal field theoryE3: they were 
able to compute the anomalous exponents for the corre- 
lation functions. Several issues remain unclear for such 
models in a magnetic field: in particular what the exci- 
tations are. Since spin-charge separation does not occur, 



the holon and the spinon cannot be the elementary exci- 
tations of the system anymore. To answer the question 
we have to turn to the low-energy effective theory. Frahm 
and Korepin's results imply that an effective description 
in terms of the gaussian model should be possible since 
conformal invariance is realized. We will find a gener- 
alization of the spin-charge separated gaussian hamilto- 
nian suitable for a description of the Hubbard model in 
a magnetic field. Our formalism is very similar to Wen's 
K matrix approach to edge states of the FQHE. This 
will enable us to characterize very precisely, in the low- 
energy limit, the properties of ID gapless models with 
or without spin-charge separation such as the Hubbard 
model in a magnetic field: we will find that in the latter 
case, although there is no spin-charge separation, there is 
still a generalized decoupling. The excitations are again 
fractional; as expected the holon and the spinon are no 
longer present in the spectrum and we will give the gen- 
eral framework allowing the description of the fractional 
states which replace them. 



A. Spin-charge separated Luttinger liquid. 

We start with the standard case when spin-charge sep- 
aration exists. Although fractional excitations are clearly 
present in Bethe Ansatz, no description of these special 
states was attempted in the low-energy limit through 
bosonization. In the following we answer several ques- 
tions: how does the holon evolve with interaction? What 
would be an effective wavefunction for it? Is it a semion? 
First, we consider the ground state of the two-component 
gaussian model, because it will suggest to us a possible 
generalization of the gaussian model which will prove to 
be the correct one for the description of gapless mod- 
els without spin-charge separation such as the Hubbard 
model in a magnetic field. 



1. Ground state of the gaussian hamiltonian. 

We consider a two component model by introducing 
an internal quantum number such as the SU{2) spin. 
We consider the charge and spin densities as well as their 
associated phase fields: 

Pc = Pi + Pi; Ps ^ PT - Pi, (3.1) 

p. = -^a,<f,, a=t,i (3.2) 

[e,(x),9,$,,(y)] ^t6„„,d{x~y), (3.3) 

'i'c/s = (3-4) 

[Qr{x),dj;^r'{y)] = iKr'^ix - y); T = C,S. (3.5) 

The effective hamiltonian derived for instance from the 
Hubbard model in the absence of a magnetic field is: 
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H = Hc + Hs, (3.6) 
Hr = ^ dx K-^ {V^rf + KriVQrf; T = C,S. (3.7) 

V$r = n,- is canonically conjugate to the field &r- One 
easily extracts the ground state which is simply a product 
of gaussians: 

T=c,s ^^-^ njiO 

(3.8) 

In terms of charge and spin densities: 

*o({PT,gJ)= n ^^P(~7^Z1 ^^l^^'«"^^'-«'•^• 



to particles with spin a) and if one sets z = expi^x, 
one easily finds: 



(3.14) 



The wavefunction is bosonic; for the fermionic LL, one 
undoes the Jordan- Wigner transformation which leads 
to: 

*F,o({a;i,(Ti}) = n l^i - ^jf 

i<j 



n 



[Zt - Zj) 
\Zi - Zj\ 



expi-sgn{ai- aj)> . (3.15) 



(3.9) 



The ground state displays of course a complete decou- 
pling of spin and charge as is apparent from the previous 
expression. This is also a Jastrow wavefunction. In real 
space: 



(The antisymmetrizing factor consists of two parts, one 
which ensures that particles of the same species anticom- 
mute, and a second part known as a Klein factor which 
allows antisymmetry for particles of different spin.) Let 
us redefine the matrix elements of g per: 



Mpr)= n e^pj^ 



4K^ 



fx fi 

K-'=X + f, 
K'=X-fi 



(3.16) 

(3.17) 
(3.18) 



// 



dxdx'pr{x) In 



sm 



it{x — x') 



pAx')]. (3.10) 



We define the charge and spin parts of the ground state 
per: 



*c/s = exp 



c/s 



dxdx' Pc/g{x) In 



sm 



7r(a; — x') 



p,,,{x')]. (3.11) 



If we denote the coordinates of particles with spin t and 
[ respectively by u and v then the ground state can be 
rewritten as (for convenience the antisymmetrizing factor 
is omitted): 

i<j i<j i,j 

(3.19) 



The previous ground state may be rewritten in terms of 
the densities of each species: 



*o({p<t}) = exp 



1, 



n 

i<j 



Ui - Ui \ \Vi - Vj \ \Ui - Vi 



1/2K, 



i<j 



(3.20) 



dxdx' pa{x)gaa' In 



. Trfx — x') 
sm 



p„,{x% (3.12) 



where we have introduced the following g matrix: 



'[\\Ur - Uj\\vi - Vj\ I \ui - Vj\ 
i<j 

\<Ti<7j/2Ke 



l/2Ks 



2 2 
2 2 



(3.13) 



(3.21) 



i<3 



The eigenvalues of that matrix are simply the inverses of 

the Luttingcr liquid parameters and K~^. If one 

rewrites the wavefunction in terms of individual electron 
coordinates pa{x) = ^ S{x — Xi) (the sum is restricted 



For the fermionic LL the charge part gets an addi- 
tional factor ^ j^^I^^ i j and the spin part, a factor 

(li^) " expifs5n(£7i-c7j). 
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These are ID Laughlin multi-component wavefunc- 
tions. In 2D they are known as Halperin wavefunctitffls 
which describe multi-component systems of the FQHEE3. 
In that context-the g matrix is known as Wen's topolog- 
ical K matrixes. The main difference between the two 
matrices is that the the entries of the K matrix are inte- 
gers while g matrix elements are arbitrary real numbers 
(only constrained to yield positive and real eigenvalues). 
The ^ matrix does not allow for a topological interpreta- 
tion either since there is no topological quantization as in 
the FQHE. We will call the 5 matrix, the charge matrix 
because it corresponds to the couplings, between particles 
in the plasma analogy (see eq.(3.12))E2l. 



2. Elementary excitations : the holon, the spinon, Laughlin 
quasiparticles. 

We generalize the approach followed for the spinless 
LL. There is a decoupling of the dynamics at two levels: 
chiral separation as well as spin-charge separation. In 
particular both the charge and the spin hamiltonians - 
He and Hg - display chiral separation: He ~ Hc+ + He- 
and Hs = Hs+ + Hg 



where the four hamiltonians all 
commute ( [_ffj,/s±: ^c/s±] = 0)- The following operators 
create the exact eigenstates of the relevant chiral hamil- 
tonians: 



V^{Qr.±Tq) = j dxexp igxexp —1^/7172 (5r,±0T, ± 

(3.22) 



6r,± = ©r T ^r/Kr] T = C,S, 

2nn 27r Ql ± 



(3.23) 
(3.24) 

(3.25) 



(The square root \/2 in the exponential comes from the 
normalization of the charge and spin fields; c and s index 
charge and spin respectively; r = ±1 for charge and spin 
respectively.) One easily checks that: 



Q^+rQ^,V^{Q-r.,±) ^6rr'Qr^±V^{Qr'^±). (3.26) 



This implies that these excitations either carry a charge 
Q = Qc,± but then have no spin (the operators VJ^), or 
that they have a spin Sz = S± — Qs.±/'2 but no charge 
(operators Fj^).As expected the fractional states come 
in two brands: the first corresponds to charge excitations 
and the second to spin excitations. Hereafter we will note 
the charge and the spins of these excitations as 



Q± Qc±, 
5± = %±. 



(3.27) 
(3.28) 



Because of the obvious relevance to physical systems we 
focus first on elementary excitations for a fermionic LL; 
we will consider the case of a bosonic LL later in section 



( III C ). Solving the constraints on the charge and current 
which again obey the selection rule — J-^ = 2n-^ and 
Ql — Ji = 2ni ( the n are integers), wc find: 





(3.29) 



This compact equation must be read as follows. Each en- 
try represents a fractional excitation; the first two lines 
are charge spinless excitations, while the last two lines 
represent spin excitations. For instance the entry 
is associated with a fractional excitation with charge 
Q — (5+ which carries no spin, and propagates in the 
right direction: therefore in the first line means 

a spinless state with charge Q = ^^2" S^i^^S 
right. Likewise the second line characterizes charge exci- 
tations propagating to the left. The line means that 
the states have no charge, a spin component S^. = 5"+ 
and propagate to the right: for instance 1/2 is a spin 
one-half fractional state. Each line gives the decomposi- 
tion of a given fractional excitation into elementary ex- 
citations: for instance the Qj^- excitation is made up of 
n| -|- ni excitations V^{Q = 1), and J| -f excitations 
V^{Q — ^^2" ) • The previous equation summarizes the 
selection rules which are obeyed by the elementary exci- 
tations. 

Let us give an example. Suppose one adds a spin up 
electron at the Fermi level in the Luttinger liquid. This 
is a Qi = 1 = J| and Qj. = = excitation or in terms 
of n| and n^, this is a (rij = 0, rij^ = 0, J| = 1, = 0) 
state. Equation (3.29) shows that the spinor 




1+ 



(3.30) 



is created; this means that the spin up electron added 
at the Fermi level kp splits into four fractional states: a 
charge ^^2° anyon propagating at velocity Uc\ this state 
has no spin; a second charge anyon with charge ^^-^^ ^"^^ 
velocity —u^ and then two spin anyons with velocities Us 
and —Us and respective spin — ■ In the special 
case of spin rotational invariance (Kg = 1) there is only 
one spin anyon: the spinon with spin ~ 1/2 which 
propagates to the right with velocity Us ( or to the left 
with velocity —Ug if the electron had been added at the 
left Fermi point —kp)- Likewise, if Kc — 1, there is 
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a single charge state, which has charge Q = 1. In the 
non-interacting case, the charge velocity Uc and the spin 
velocity Ug are equal and therefore the spinon and the 
charge 1 state subsume into a single state since they move 
in the same direction with the same velocity: we have just 
recovered the spin up electron. 

a. Holon and spinon. Let us identify the content of 
the elementary excitations, starting with V^{Q = 1) 
and V'^{Sz = 1/2). The charge and the spin carried 
by these fractional excitations make it reasonable to in- 
terpret them as the holon and the spinon respectively. 

The physical processe involved in the creation of each 
of these states confirms this identification. Indeed the 
minimal operation which involves V^{Q = 1) is obtained 
when J| = = and set n-f = 1 = in equation ( 3.29 ). 
This is an excitation for which — I — Qi and Jf = 
= J| which means that this is a pure charge process (no 
spin variation Sz = (Qt ~ Ql)/'^ = Oi spin current nor 
charge current). V^{Q — 1) is an excitation associated 
with the addition of charge in the LL. All the transitions 
in Fock space which occur after adding a charge to the 
ground state therefore involve V^{Q = 1)- rJ^ is then 
consistent to identify VJ^((5 = 1) as the holonE3. 

Likewise the minimal excitation generating V^{Sz = 
1/2) is a spin one transition which is a pure spin process. 
All excitations for which there is a spin flip will therefore 
create V^{Sz — 1/2) (in pairs). This is what we expect 
from a spinon. 

Notice that both for the holon and spinon there are 
even-odd effects arising in the low-energy gaussian the- 
ory. Indeed, equation ( [3.29| ) shows that an excitation 
with includes a spin one-half transition will not create a 
spinon (we have instead a "hybrid" spin excitation): one 
needs at least a spin one transition to crate a spinon. 
This means that adding a single electron does not create 
a spinon: even number of electrons are required. This 
makes sense since for a spin chain the minimal spin ex- 
citation is also a spin-flip which involves two electrons 
and not just one. The same behaviour is observed with 
the holon: the minimal process which creates it adds two 
electrons {Q = 2 since Qi — 1 — Qi). These even-odd 
effects are a direct consequence of statistics and would 
not be observed with a bosonic two component LL. 

With eq.(pll), we find that the wavefunctions for the 
holon and the spinon are simply (with z — expi^x): 



'^hoionizo) - V+{Qc = 1, Zo)*c 

1 



n 



ni 



l/2K^ 



n 

i< 

n( 



i<j 



n 



l/2Ka 



n 

l<j 



xo.i<j 

Zi ~ Zj\ 



(Xj — Xj) 
\Xi — Xj I 



1/2 



X exp - 



2kF 
I 



N 



Xq 



n 



[Xi Xj) I 



, (3.32) 



and: 



^spinon (cq, Zq) 

= V+{Sz^<Jo/2,zo)-^s 



n 



\xi - xo\ 



ni 

l<j 



n 

i 



\Zi - Zo 



cjQaij2Ks 



xo,i<j 



{•^i -^j ) 
\xi — Xj I 



CTiO-j /2 



(3.33) 



The holon and the spinon are both anyons with ex- 
change statistics: 



TT 

2i^' 

TT 



2K. 



(3.34) 
(3.35) 



(the statistics were computed in subsection IB). Except 
for the special case Kr — 1, these objects are not semions; 
in addition, for holons we must also require Vc ^ Vg to 
ensur e spin-charge separation. Contrast our results (eqns 
( |3.31 ) and ( 3.33| )) with the commonly used but ipcorrect 
charge-spin decoupling of the electron operator^. The 
holon and the spinon generalize the dual excitation found 
for the spinless LL; in the same way, duality transforms 
the holon and the spinon into the two-component gener- 
alizations of the Laughlin quasiparticles. This is a most 
remarkable yet simple result because it shows that two 
seemingly unrelated fractional excitations -the holon (or 
spinon) and the Laughlin quasiparticle- occurring in two 
very different contexts are actually deeply connected. 

Similarly to the spinless LL, in addition to the holon 
and to the spinon, we hybrid excitations complete the 
basis of fractional excitations. Their charge and spin are 
intermediate between those of the holon and spinon and 
those of their dual excitations, the Laughlin quasiparti- 
cles which we discuss now. 

b. Laughlin quasiparticles. We can choose another 
basis of elementary excitations dual to the previous one 
which will parametrize the excitations in terms of cur- 
rent processes and electron addition at the Fermi sur- 
face. This basis, emphasizing Laughlin quasiparticles as 
elementary excitations reads: 



(3.31) 



/Q+\ 




Q- 
s+ 




\S- J 






(3.36) 
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J 



2n-\ and Qi — Ji = 2n|. In 



where again Qi 

addition to the hybrid quasiparticles V^{Q — ^^^) 
and V^{Sz = ^^^) which already existed in the pre- 
vious basis, we have two excitations associated with pure 
charge current or spin current processes: Vc{Qc = Kc) 
and Vs{Sz = Ks/2). Under electromagnetic duality the 
latter are conjugate to the holon and spinon respectively. 
Actually they are obtained by spin-charge separation of 
the two-component Laughlin quasiparticle, and we may 
call them a Laughlin holon and a Laughlin spinon. Let 
us compute the wavefunctions of these two excitations; 
one finds: 



i 

X exp -i2kF + ^0 ] *c, 



(3.37) 



i 

X exp-i2ao(fcT~^i)f^^^I^^+a;o)«'s (3.38) 



M 



(3.39) 



In the previous expression do takes on the values ±1 
and fc| , ki are the Fermi vectors associated with particles 
of spin up and down: fco- = J-^Vg- and N — N-^ — N^, 
M = — N j . 'I'c and ^'s are given for the bosonic 
LL in eq.( ^.2C| ) and ( |3.21 ); for the fermionic LL, there 
are additional phase factors given in the text follow- 
ing eq.( [3.20 ) and ( [3.21 ). Plasma analogy allo ws t o get 
the charge and spin of the Laughlin holon eq.(3.37) and 
spinon eq.(3.39)). The Laughlin holon and spinon are 



the stable excitations into which the Laughlin quasipar- 
ticle decays as a result of spin-charge separation; in the 
localized Wannier basis we have considered throughout, 
the product of the wavefunctions of the two excitations 
yields indeed: 



,(Z0, (To) = Y[{Zi - ZoY'O'^ Yi 



(3.40) 



which is just the generalization of the Laughlin quasi- 
particle to two component systems: when we add spin, 
the Laughlin quasiparticle comes in two flavours (up or 
down) and the Laughlin correlation hole acts only on 
particles of the same flavour. It therefore carries both 
fractional charge and fractional spin. (It would be spin- 
less if the Laughlin prefactor were J^,j(zi — ^o) instead of 

n,(z,-2o)'"'';0- 

The Laughlin holon and spinon have statistical phases 
9c — ttKc and 6s = TrKs/2. For a spin-rotational invari- 
ant system, the Laughlin spinon and its dual conjugate - 
the spinon - are identical states {Kg = 1 is the self-dual 
point for spin excitations). 



B. Luttinger liquid without spin-charge separation. 

1. The general LL and the charge matrix. 

We now generalize the standard LL theory to include 
situations with no spin-charge separation. We start from 
the ground state of the gaussian hamiltonian: 



*o [9] = exp 



dxdx' Pa{x)gaa' In 



n(x — x') 
sm 



(3.41) 



(3.42) 

(3.43) 
(3.44) 



and relax the constraint — (while g is kept sym- 

A ^J, 



metric) . We consider the charge matrix: g^^' = 



^ A' 



and the associated wavefunction 5*0 [5] . We introduce for 
convenience the eigenvalues of the charge matrix and the 
unitary matrix P : 



P^^gP = D = 

5t7 



1/Ki 







I/K2 



(3.45) 



Kr 



(3.46) 

The normal modes of the charge matrix are simply: 

Pr = Par Pa <^ Pa = ParPr T = 1, 2 (3.47) 

^Pa{x)gaa'Pa'{x') = ^ Pt{x) Pr{x') . (3.48) 
aa' T ^ 

If the charge matrix obeys a Z2 symmetry then the nor- 
mal modes are just the charge and spin density (up to a 
normalization factor) pi = pc/V2 and p2 = Ps/V2). The 
wavefunction now reads: 



exp 



2Ki 



X exp 



2Ko 



dxdx'pi (x) In 



dxdx' p2{x) In 



sm 



'k{x — x') 



nix — x') 
sm 



Pi{x') 



P2{X') 



(3.49) 



We have expressed the ground state in this decoupled 
form because this allows us to directly write down a gaus- 
sian hamiltonian with ground state ^'0(5]. This general- 
izes the spin-charge decoupled gaussian theory. We intro- 
duce the phase fields associated with the normal densities 
pi and p2 
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--^dx^r, jr = A=dxQT, 



[^r{x), d^Qr'iy)] = i5rr'S{x - y). 



(3.50) 



(3.51) 



It is then clear that "io \g\ is the exact ground state of the 
following family of two component gaussian hamiltonians 
for arbitrary velocities ui,U2'- 

H\g,ui,U2] 

= Hb[ui,Ki]+Hb[u2,K2] 



r=l,2 ^ -^0 



dx 



K-^{d,^rf + Kr{d,Qrf ■ (3.52) 



The two hamiltonians HbIui^Ki] and Hb[u2, K2] com- 
mute by construction. The next section will be devoted 
to the properties of that generalized LL theory. Let us 
stress here the main property of this general LL now: 
by construction that theory corresponds to a generalized 
separation: the normal modes will not be charge and spin 
modes but mix charge and spin in a proportion fixed in 
time. This will translate for the fractional excitations to 
states with both fractional charge and fractional spin. 



Anomalous exponents are easily computed as functions 
of the charge matrix 5 which leads to compact expressions 
valid both for the LL with spin-charge separation or for 
the more general LL; one introduces the Fermi vectors: 
kpa = ^^jf-- For instance density-density correlators are 
obtained with H[g] and with the bosonization formulas, 
yielding the static structure factor; the dominant Fourier 
components are k = 0, k = 2kF-[, k = ^kpi and k = 
2fcFT + 2fcFi. Near fc = 0: 



2{nx)^ ' 



< Sp{0)6p{x) >-- 



{ttx)^ ' 



fe=0 



E 



9 ^aa' 



(3.57) 



(3.58) 



For the higher harmonics one includes a mode at 
2kFi + '2kFi = 27rp (which appears in the Hubbard model 
in a magnetic field): 



(7r.T)2 

cos{2kF^ + 2kFi)x 



< 5p{0)Sp{x) > 
+b 



cos{2kF-ix) cos2fci?4a; 



+ ai 



X 



2+a(2kFT+2kFi) 



j.2+a(2kFi) 
(3.59) 



2. Main properties. 

The compressibility and spin susceptibility are easily 
computed and one finds: 



1 d^Eo TT 



_i ^ 1 d^Eo 



(3.53) 



(3.54) 



po and ps are the charge and spin mean densities, while 
Eq is the ground state energy. The Drude peak is: 



D=- 



1 d^Ep 



(3.55) 



where ^ is a flux threading the LL ring of length L. 
This expression can also be recovered using the Kubo for- 
mula; one then needs the expression of the current den- 
sity which one finds with the continuity equation. The 
current is renormalized as for the spinless LL and the LL 
with spin-charge separation. In contrast to the case of 
the LL with spin-charge separation, here the expression 
involves both Ki and K2 because both modes one and 
two involve charge: 



jR{x) 



E 



(3.56) 



2 + a(y2kFa) 
a(2kF\ + 2kFi) 



(3.60) 
(3.61) 



Ak=o is fixed in the low-energy limit but the other 
constants a-f ,a| , h are non-universal and depend on high- 
energy processes. When there is spin-charge separation, 
the exponent for AkF oscillations and the constant Ak=o 
are related by the equation 2 -|- a{AkF) = AAk=o; in the 
general case, we find: 

2 + a(2fcFT + 2kFi) = 4Afc=o 

The derivation of the exponents is done in exactly the 
same manner as in the spin-charge separated LL. 
Spin-spin correlation functions are: 



< SMSz{x) >= 

^ COS 2kaX 



J2aa' (o-O-'ff Ka') 



2{Trx) 



\x 



< S^O)S-{x) >= '^^(kF, +kFi)x 

\x\ 

9n-9n + I^i9~^ +9)a 



7 = 



Electronic Green functions decay as: 



(3.62) 

(3.63) 
(3.64) 
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\x\ 

(3.65) 

For bosons the exponent is modified as: 1 + as(cr) = 
^gaa- These exponents are derived with the bosonization 
formulas but can also be found by plasma analogy. 



3. The charge matrix: a summary. 

For the two-component LL there are three interesting 
situations which we summarize below: 

i) in the general case, the (symmetric) g matrix has 
arbitrary entries; there is no spin-charge separation but 
a more general two-mode separation: 



9 



A /i 
A' 



(3.66) 



As will be shown below, the Hubbard model in a mag- 
netic field can be described by such a theory. 

ii) the g matrix has a Z2 symmetry; this case pertains 
to spin-charge separation: 



1 



A — /i 



(3.67) 



Indeed the symmetry under the exchange of up and down 
spins implies that the normal modes of the charge matrix 
are just the charge and spin modes. The LL parameters 
are then the eigenvalues of the inverse of the charge ma- 
trix. This situation describes models with spin-charge 
separation but with a spin anisotropy, for instance a 
Hubbard model to which one would add some Ising term 

iii) the g matrix corresponds to a SU (2) symmetric 
case {K„ = 1): 



9 



^ + 1 



1 



2fi + l 



1 (3.68) 



This situation describes the low-energy limit of the Hub- 
bard model. It is noteworthy that the wavefunctions 5" \g\ 
for that sub-case were used in a variational approach of 
the \T) t — J model giving very good results although 
it was not realized tiijey were the exact ground states of 
the gaussian modelo. The reason why it is so is now 
transparent. 



C. Elementary excitations for the generalized LL. 

We now consider the excitations of the general bosonic 
and fermionic LL with or without spin-charge separation. 
Let us inject particles in the Luttinger liquid. In real 
space this is described by the operator: 



V{x) = exp -iv^ V {QaQa[x) - J,$,(x)) . (3.69) 



Fractionalization stems from two decouplings: chiral sep- 
aration and a separation for the internal quantum num- 
ber generalizing spin-charge separation. In terms of the 
normal modes fields Qr and <I>r (r = 1, 2): 



V = 



exp-zV^^ I (^p,,g^)er-(^p,.j,)$r] 

r \ (7 a J 



The chiral fields are: 

Qr,±{x) =Qr(x)T'^r(x)lKr, 

and therefore: 

y{x) = JJexp-i\/7''Qr,±0r,±(a;)- 



(3.70) 
(3.71) 
(3.72) 



This expression explicitly shows a decoupling into four 
components. We have defined in the above the chiral 
charges: 



(^ Pa.Qa) ± ifr(^ Par J.) 



(3.73) 



The following operators are exact eigenstates of each chi- 
ral hamiltonian -ff±,r t = 1,2: 

K-^(Qt,±, 9) = y dxexpiqxexp— «A/7r<3r,±6r,±(a;), (3.74) 
27m 27r Qr.± 



9 = 



L L Kr 



(3.75) 



The chiral charges correspond to the charge and spin car- 
ried by each of these excitations up to a normalization 
factor: 



..V^{Qr.,±.q) \ ^Qr.±ParV^{Qr.±,q), (3.76) 



which implies that the charge and spin of V^{Qr,±Tq) 
are: 



(3.77) 
(3.78) 



Thus for an arbitrary charge matrix, fractional excita- 
tions carry both charge and spin. However the ratio of 
charge to spin is constant for each given mode r = 1,2: 
Scr <^^crr) Q = 'S'z (X^cr -Pcrr); of coursc the phouous 
associated to each mode mix charge and spin in exactly 
the same proportions since: 
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Prix) = (^^^P-yrj Pc{x) + {]^^(^P<yr^ Ps{x) 
= (^E^-) Pc{x)+{^P.r 



s^{x) (3.79) 



where Sz{x) is a spin density. 

It is convenient to define the charge to spin ratio: 



Q_ 

2S' 



(3.80) 



for each mode. Unitary implies that if for the first mode: 

Q 



2S, 

then for the second mode: 



P 



p' 



(3.81) 



(3.82) 



where p is arbitrary. Note that for a Fermi Hquid these 
ratios are r = ±1 (we are characterizing Landau quasi- 
particles (or holes) of either spin) and when spin-charge 
separation is realized the ratio is either r = or r = ±oo. 

Let us give the elementary excitations. The simplest 
case is that of bosons: 



Qt,± Qt 



Ql 



Plr 



(3.83) 



To simplify the notation, we have only written a sin- 
gle line, but Qt,± and the other entries should be read 
as four-vectors. Q„ and ^ are arbitrary independent 
integers, which shows that the states V^{QT,±Tq) where 
Qt,± = '^P^tKt or ztPj^^i^T- are elementary ex- 

citations. As an illustration let us consider the simple 
case of a Z2 symmetric charge matrix for bosons which 
have a pseudo-spin index. The unitary matrix P is: 



P. 



l/x/2 1/V2 \ 

1/V2 -i/V2y 



(3.84) 



Then it follows from eq.( 3.77D and ( [3.78 ) that for mode 
T = 1 (the charge mode) Q — \/2(3r.± and Sz — for 
mode T = 2 Q = but Sz — Qt,±/ \/2- If we take these 
normalizations into account: 

(Q+' 
Q- 




Once again we find a charge 1/2 particle and a charge 
Kc Laughlin quasiparticle as for the bosonic spinless 
LL. But in addition we find new states resulting from 
a fractionalization of " pseudo-spin" for the bosons (half- 
spinons for instance). 

For the fermionic LL the elementary excitations are 
obtained by the equation: 



l±Kr\ ( l±K. 



2 J • y ^' 2 
+n^{TP]rKr) + n^{TPirKr). (3.86) 

where we have resolved the constraint: Qa — — ^n^. 
This fully characterizes the low-energy elementary exci- 
tations of a LL in a magnetic field (see below). 



D. Application to the Hubbard model. 

To illustrate the previous results we discuss the Hub- 
bard model in one dimension. The model was solved ex- 
actly by Bethe Ansatz by Lieb and Wu. In zero magnetic 
field, for repulsive {U > 0) interactions, a LL metallic 
phase exists both for weak and strong coupling, except 
at half-filling. For very large U the spin-charge decou- 
pling is valid at all energy scales. This was shown by 
Ogata and Shiba who also found that the Bethe Ansatz 
ground state then took a remarkable factorized formE3: 
it is the product of a charge part (a Slater determinant 
for free fermions involving all electrons ) and a Bethe 
wavefunction similar to that of the Heisenberg model on 
a reduced lattice from which one has removed the holes 

"^HubbardiXi^CTi) 
= det{e^p ikjn, \kj\ < kp) Heisenberg{Vi,'7i) (3.87) 

{yi is the coordinate in the reduced lattice of particle i 
whose real position is Xi). 

It is instructive to compare it to the two-component 
Jastrow wavefunctions which are also explicitly spin- 
charge decoupled. The Slater determinant is rewritten 
as (in terms of the circular coordinates z): 



'^Hubbard{Xi,(Ji) — JJ^ (^i ^ ^j) HeisenbergiVi, C^i) 

(3.88) 



This is to be compared with: 



I 



<5<7,<T 



n \ ir^^) exp^;.gn(a, - a,) } (3.89) 



(3.85) Or if we separate spin and charge: 
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/2Ka ( i^i - Zj) 



1/2 



(3.90) 



i<j 



— ^ expi-s5n(cri - cTj)]. (3.91) 

The spin part of the Laughhn ground state is just the 
Haldane-Shastry wavefunction if = 1 (rotational in- 
variance) which has the same large-distance physics as 
the Heisenberg ground state. We can also determine Kc 
without any computation by just reading off its value 
from the wavefunctions: the charge parts of the two wave- 
functions coincide if Kc = 1/2 which indeed is the known 
value of the LL parameter for large U. 

Bethe Ansatz gives the spectrum and the eigenstates; 
however it is very difficult to compute correlation func- 
tions. An important advance came however with the 
works of Frahm and Korepin who used CFT in conjuiic.- 
tion with Bethe Ansatz to compute critical exponentsc3. 
If a theory is conformally invariant, one can show that the 
finite-size energies of excitations are directly related to 
their operator dimension (which is one-half of the anoma- 
lous dimension of their correlation function). By using 
Woynarovich's Bethe Ansatz calculations for the finite- 
size spectrum to order 1/L which he computed within 
a so-called "dressed charge matrix formalism"LJ, Frahm 
and Korepin were able to extract critical exponents for 
the correlation functions of the Hubbard model. In par- 
ticular they found that in the presence of a magnetic 
field, spin-charge separation was not realized. Penc and 
Solyom later showed that in 1/L the same spectrum de- 
rived by Woynarovich could be expressed in terms of 
a generalized Tomonaga-Luttinger model with interac- 
tions described in the g-ology framework; using equa- 
tion of motion methods they also derived the anomalous 
exponentsEJ. These two approaches give little insight into 
the nature of the elementary excitations: how are the 
holon and spinon modified as a function of microscopic 
parameters? The description of spin-charge separation 
(or its absence) is not transparent either: the dressed 
charge matrix tell us little about spin-charge separation; 
its changes are not easy to relate to that property. This 
is to be contrasted with our charge matrix formalism in 
which spin-charge separation is directly connected to a 
symmetry of the charge matrix g {Z2 symmetry) . We will 
show that the " dressed charge matrix" of Bethe Ansatz 
and the charge matrix g are in fact related: the inverse 
of the symmetric charge matrix is roughly the square of 
the Z matrix. We will proceed in the following manner: 
we will show that Woynarovich's finite-size spectrum is 
identical to that of our generalized LL. This yields the 
charge matrix g in terms of the dressed charge matrix Z 
and gives us both the anomalous exponents and the frac- 
tional excitations since we already derived them for the 
generalized LL. The relation of our charge matrix formal- 



ism to Penc and Solyom g-ology approach is the follow- 
ing: it can be understood as a bosonization of their gen- 
eralized Tomonaga-Luttinger model; it is much simpler 
however to work directly within the gaussian hamiltonian 
framework. Our approach has several advantages in ad- 
dition to making an explicit contact with the seemingly 
unrelated physics of Laughlin states: (a) we avoid an am- 
biguity in the determination of aaomalous exponents in 
Frahm's and Korepin approach;t3 (b), we can give the 
nature of elementary excitations ( and show that they 
are fractional states in the first place ) and (c) we are 
able to give a clear criterion of spin-charge separation. 

Woynarovich's finite-size spectrum in Frahm and Ko- 
repin's notations is the following: 

E{^N,D,Nt,Nt)-Eo 

= ^ [vc{Nt + N-) + Vs{N+ + N-)] 



27r 



-AN'^{Z^^fVZ-'^AN + D'^ZVZ'^D 



P{AN,D,Nt,Nt)^Po 

V(iV+ -N-) + Vs{N+ - N~)] 



(3.92) 



L 



27r 



[AN^D] + 2DckF] + 2{Dc + Ds)kFi. (3.93) 



kp^ — and kpi = are the Fermi momen- 

tum for particles of spin up and spin down. The energy 
and the momentum are those of a state with the (inte- 
ger) quantum numbers {AN, D,N^ , N^); there are two 
modes indexed by c and s: these two modes do not in 
general correspond to charge and spin. Z is a 2 by 2 
matrix: 



Z = 



Zee Ze. 
Zsc Zs. 



(3.94) 



and AiV and D are two-vectors;-, AA^ = {Ne = + 
A^X, A^ = Ni) and D = {De,Ds)M In these expressions 
are integers: they are simply the modulus of phonon 
momenta in units of 2tt/L for the two modes c and s; the 
index ± refers to the sign of the momentum. The phonon 
velocities for the two modes are (vc, Vs). 

The spectrum of the general gaussian model H [ur,g] 

is: 

= ^ K = l(iV,ti + A^-.i) + Vr=2{N+^2 + K=2}] 



01 



(3.95) 



p(g,, j,,A^±) 
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ETTt 



(3.96) 



iV^j^,iV^2 again the moduli of phonon momenta. 
The charges and currents {Qt, Jt) are related to {Q^, J a) 
by Qr = P^T-Q^ and Jr = Par J a- We can now identify 
the parameters of both theories: 



Vr = Vc/s- 



N 



(3.97) 
(3.98) 
(3.99) 



The zero modes can be identified term by term; it is 
sufficient to consider the current terms to uniquely de- 
termine the charge matrix. The charge zero modes yield 
extra relations which lead to the very same expression 
for g. Indeed expanding the squares gives; 



Since: 



KiPf^ = 

KiPfi = 
KiPnPli = 
K2PI2 = 

K2P^2Pl2 = 



9 



{Zee ~ Zse)'^, 
{Zse)\ 

{^cc ^ sc)Z SC-i 
{Zss)\ 

{Zcs ~ Zss)Zss- 



K-j-Pfj-YPa' r 



(3.100) 



(3.101) 



it follows that the inverse of the charge matrix is: 

9^^^] = {Zee ^ Zsci^ + {Zes ^ Zssf' 

= {Z^ef + {Zssf 
9 ^ 9 "'"XI ^ {Zee ^ Zsc)Zsc + {Zes ^ Zss)Zss (3.102) 

We define the matrix Z obtained from the dressed charge 
matrix Z by subtracting the second line from the first: 



Z = 



Then: 
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Z sc Zes Zg 
Zss 



zz' 



(3.103) 



(3.104) 



This is the most important result of the present sec- 
tion: the low-energy properties of the Hubbard model 
are expressed in terms of quantities which can be com- 
puted from the microscopic parameters and spin-charge 
separation simply follows from the Z^ symmetry of our 
charge matrix. In the framework of the dressed charge 
matrix Z approach, the second feature is not easily de- 
coded from the structure of Z which is then triangular 
with some relations between its matri&Telements whose 
physical interpretation is quite unclear .E3 We can use the 



results of the previous sections on the elementary exci- 
tations and those on the various properties of the charge 
matrix hamiltonian such as the Drude peak, the suscep- 
tibility, the anomalous exponents. In particular the new 
modes replacing the spin and charge modes are simply 
the eigenvectors of the charge matrix. 

The charge matrix is obtained by inversion: 



1 



9 



9 \i 



(detZ)' 



-9 \i 



9 n -9 TT 



(3.105) 



Term by term identification of the charge zero modes 
Qa would lead to exactly the same expression for in- 
deed: 



^^^^^^^ 



{Zss? 

(dot Zf 

{Zes — 



(det zy 

^{Zes ^ Zgg)Zs 



(det zy 



{Zsef 



and since 



P\2Pi2 - 



(detZ)^ ' 

{Zee ^ Zge? 

(det zf ' 

^{Zee ^ Zse)Zg. 

(det zf 



(3.106) 



. _ _r' = Y^T^T^ParPa'T One rccovcrs 
equ. ( |3.105D . As it should be, one can check that the 
anomalous exponents predicted for H [ur , g\ agree then 
completely with Frahm and Korepin's results. 

Let us illustrate these results in two situations, one 
with spin-charge separation, the other without. From 
these, we can exhibit the criterion for spin-charge sepa- 
ration within the matrix formalism. 

In the presence of spin-charge separation, Frahm and 
Korepin find that the dressed charge matrix Z va Z — 



i Ze 

^ e/2 Zg 





1/V2 



inverse g ^ are: 



9 ■ 




This implies that 5 and its 

(3.107) 
(3.108) 



The charge matrix explicitly exhibits spin-charge separa- 
tion and takes the form characteristic of SU{2) symme- 

try. The eigenvalues of g are Ke 



= ^ and Kg = 1. 
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The Z matrix can also be explicitly computed in the 
limit of infinite repulsion with a magnetic field close to 
the critical field he for which all the spins are polarized 
(i.e. close to the ferromagnetic phase). In terms of the 
parameter 



(5 = 



he 



the dressed charge matrix is: 

^ - I 2^ ^ _ 1 5 



(3.109) 



(3.110) 



which implies that the (inverse of the) charge matrix g 
is: 



9 



(3.111) 

This expression shows explicitly the breakdown of spin- 
charge separation (except for 5 = 7r/4). 



IV. CONCLUSIONS AND PERSPECTIVES. 

The main goal of our paper was to establish and 
describe fractional excitations for the Luttinger liquid 
within the bosonization scheme, bf Bethe Ansatz gives 
exact eigenstates and shows the existence of some frac- 
tional excitations: however their description is quite com- 
plex in that framework and it is unclear how to generate 
systematically a complete set of excitations. In the low- 
energy limit, the Luttinger liquid approach allows a very 
precise characterization of the fractional states already 
known from exact solutions but what's more allows us 
to discover novel fractional excitations which may carry 
irrational charges (the ID Laughlin quasiparticle, the hy- 
brid state). In section || the low-energy spectrum of Lut- 
tinger liquids can be reinterpreted in terms of fractional 
states: for instance, the particle-hole continuum consists 
of a Laughlin quasiparticles-quasihole continuum. The 
quasiparticle perspective clarifies many properties of the 
LL: the renormalization of the current operator is a di- 
rect consequence of fractionalization; for spin chains, we 
present the correct description of the spinon excitation in 
the generic case of a violation of SU (2) invariance. We 
also show that the Sz ^ Q continuum of spin chains in- 
volves the analogs of Laughlin quasiparticles. In section 
III we describe fractional excitations such as the holon 
or the spinon for the Luttinger liquid with spin; we also 
present a generalization of the gaussian theory valid for 
Luttinger liquids without spin-charge separation and dis- 
play in that situation the new f ractional states replacing 
the holon and spinon (see |III q ). 



An important test, of course, would be to observe ex- 
perimentally (or numerically) all these fractional states. 
Although the existence of the holon and the spmon was 
ascertained theoretically quite a long time agoQ no ex- 
periment has yet allowed their detection: in fact, the 
property of spin-charge separation itself is not yet es- 
tablished experimentally. The observation of two of the 
fractional states discussed in this paper would be partic- 
ularly important: the LL Laughlin quasiparticle and the 
hybrid state. Indeed they may assume irrational charges. 
The precise spectroscopy of fractional excitations we have 
done in this paper allows us to determine which processes 
are involved in their creation: for the Laughlin states, 
current probes are needed, while the hybrid particle is 
created by addition of an odd number of electrons. For 
Laughlin quasiparticles shot noise is likely an adequate 
probe: the shot noise coefficient for Luttinger liquids cae. 
be computed exactly and is predicted to be equal to K f3 
in the two-dimensional electron gas at filling j/ = 1/3 this 
yields a charge 1/30. The latter situation involves Wen's 
chiral Luttinger liquid. The identification of the shot 
noise coefficient with the charge of a carrier has been de- 
bated because the coefficient one measures might actually 
be the conductance rather than a quasiparticle charge (at 
J/ = 1/3 the conductance also assumes the value 1/3) . For 
the non-chiral Luttinger liquid our spectroscopy of frac- 
tional states allows to resolve that ambiguity: the shot 
noise coefficient is indeed identical to the conductance K 
of the LL but the (backscattering) current-current cor- 
relation function measured in shot noise involves charge 
K excitations, because charge K LL Laughlin quasiparti- 
cles are precisely generated by current excitations. These 
might thus be detected in any physical realization of the 
Luttinger liquid: quantum wires, or possibly nanotubes. 
An intriguing possibility is also suggested by recent ex- 
periments of tunneling-at the edge of a two-dimensional 
gas in a magnetic fieldEil. I —V characteristics measured 
at the edge showed very surprising non-Fermi liquid be- 
haviour compatible with a chiral LL with unquantized 
LL parameter K] the I — V curves evolve smoothly when 
one varies the filling fraction and do not show a plateau 
structure. There seems to be a continuum of Luttinger 
liquids living at the edge: this caused quite a stir be- 
cause the chiral LL theory can presumably be derived 
only for incompressible filling fractions. These puzzling 
results are sofar unexplained, but if an interpretation in 
terms of a single-boson mode chiral LL with unquantized 
parameter K can be in some manner justified, according 
to the results given in the present paper this would im- 
ply that there exists charge K Laughlin quasiparticles in 
that experimental setting: such a chiral LL is identical 
to the chiral half of the non-chiral gaussian hamiltonian 
considered throughout our paper. 
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APPENDIX A: DISPERSION OF THE 
FRACTIONAL STATES. 

We show here that Vg^(q)|*o > (where l^^o > is the 
interacting ground state) is an exact eigenstate of the 
chiral hamihonian H± with energy: 



E{Q±,qn) = 



u \qn\ + 



TTU > 



2L K 



(Al) 



Let us rewrite the state considered VJ^^(q„)|^o > in 
t erms o f zero modes and phonon operators (we use equ. 
( ITTel) , (p7| ), ([l26|) and (|l27|)): 



eigenstates with the same eigenvalue of the appropriate 
chiral hamiltonian (_ff_|_ or This is enough to prove 

that VQ^{qn)\^Q > is an exact eigenstate of H±: 



J2u\q\a^aq + 

Lq>0 



TTU Qj. 



(A5) 
(A6) 



27m 



L 



TTU 



K?±(9")l*o>, (A7) 



where g„ 



2^ is the momentum due to phonons. 



Vq^ (2^)1*0 >= expiQi 



exp-i\/7r(3± ^ 6±^„expi— ^a;|*o > (A2) 



Taking into account the fact that the operators bq an- 
nihilate the ground state, it follows that: 



F+^(a;)|*o >-expiQ4 



-—Q+x-V^iQo-^ 



exp-iV7rQ+ ^ J^;^^6+expi^a;|*o >, (A3) 



n>0 



Ktt \n\ 



with a similar expression for Vq (a;)|5'o > (the sum is 
then over negative momentum phonons). Going back to 
reciprocal space: 



1 f Stt 

^Q±(*i)|*o > = ^ / dxcxp-i—nx 

X exp iQ 



exp-i\/7r(5+ 




^ 6+expi^a;|*o >, (A4) 



K-K \p\ 



which shows that Vq^qn)]^!) > only spans chiral 
phonons with momentum ±n > ; in other words this 
state is obtained by the action of the zero mode of the 
chiral field plus the creation of phonons with momenta of 
the same sign. When one expands the phonon exponen- 
tial, the integral over position will select configuration of 
phonons with identical total momentum q^ — AH 
these configurations consist of phonons of identical chi- 
rality and total momentum which means that they are 



1 R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983) 

^ L. Saminadayar, D. C. Glattli, Y. Jin, B. Etienne, Phys. 

Rev. Lett. 79 2526(1997); R. de Picciotto, M. Reznikov, 

M. Heiblum, V. Umansky, G. Bunin, D. Mahalu, Nature 

389, 162 (1997) 
^ J. Hubbard, Proc. Roy Soc. A 240, 539 (1957); 243, 336 

(1958); 276, 238 (1963). 
*H. A. Bethe, Z. Phys. 71, 205 (1931). 

^ E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 (1968). 
•^L. D. Fadeev, L. A. Takhtajan, Phys. Lett. A 85, 375 
(1981) 

J. D. Johnson, S. Krinsky, B. M. McCoy, Phys.Rev. A 8, 
2526 (1972) 

® F. D. M. Haldane, Phys. Rev. Lett. 45, 1358 (1980); ibid, 

J. Phys. C 14, 2585 (1981) 
^L. D. Landau, Sov. Phys. JETP 3, 920 (1957); 5, 101 
(1957); 8, 70 (1959). 

^° J. M. Luttinger, J. Math. Phys. 4, 1154 (1963). 

" S. Tomonaga, Prog. Theor. Phys. 5, 544 (1950). 

D. C. Mattis and E. H. Lieb, J. Math. Phys. 6, 304 (1963). 

" C. Bourbonnais, L.G. Caron, Int. J. Mod.Phys.B 5, 1033 
(1991); R.Shankar, Rev. Mod. Phys. 66, 129 (1994) 
We call these excitations describing collective density fluc- 
tuations interchangeably phonons, or plasmons. 
c is a real number known as the conformal anomaly. 
A. A. Belavin, A. M. Polyakov, and A. B. Zamolod- 
chikov, Nucl. Phys. B 241, 333 (1984); C. Itzykson and 
J.-M. Drouffe, Statistical Field Theory, Cambridge Uni- 
verstity Press, Cambridge, 1989, vol. 2; P. Ginsparg, in 
Fields, Strings and Critical Phenomena, Les Houches Sum- 
mer School 1988, Eds. E.Brezin, J. Zinn- Justin (North- 
Holland, Amsterdam, 1990). 

H. W. J. Blote, J. L. Cardy and M. P. Nightingale, Phys. 
Rev. Lett. 56, 742 (1986) 

N. Kawakami, S. K. Yang, Phys.Rev.Lett. 67, 2493 (1991) 
P. Bouwknegt, A. Ludwig, K. Scoutens, Phys. Lett. B 338, 
448 (1994); K. Schoutens, Phys. Rev. Lett. 79, 2608 (1997) 
2" O. Babelon, H. de Vega, and CM Viallet, Nucl. Phys. B 
220, 13 (1983); K. Sogo, Phys. Lett. A 104, 51 (1984). 
F. D. M.Haldane, Phys. Rev. Lett. 60, 635 (1988); B. S. 
Shastry Phys. Rev. Lett. 60, 639 (1988) 



25 



F. D. M. Haldane, in Proceedings of the 16th Taniguchi 
Symposium on Condensed Matter Physics, Kagoshima, 
Japan, 1993, Eds A.Ojiki, N.Kawakami (Springer, Berlin- 
Heidelberg-New- York, 1994) 

^■^ We omitted a phase factor corresponding to a momentum 
n boost. It can be ignored after a bipartite transformation. 
F. Woynarowich, J. Phys. A 15, 2985 (1982) 

25 H. Frahm, V. E. Korepin, Phys. Rev. B 42, 10553 (1990); 
ibid, Phys. Rev. B 43, 5653 (1991) 

F. Calogero, J. Math. Phys. 10, 2191(1969); B. Sutherland, 
Phys. Rev. A 4, 2019 (1971) 

M. R. Zirnbauer, F. D. M. Haldane, Phys. Rev. B 52, 8729 
(1995); 

2**X. G. Wen, Phys. Rev. B 41, 12838 (1990); X. G. Wen, 

Adv. Phys. 44, 405 (1995) 
2^ Z. N. C. Ha, Phys. Rev. Lett 73, 1574 (1994) 

M. P. A. Fisher, L. I. Glazman, in Mesoscopic electron 

transport, ed. by L. Kovenhoven, G. Schoen, L. Sohn 

(Kluwer, Dordrecht). 

R. Heidenreich, R. Seller, and A. Uhlenbrock, J. Stat. Phys. 
22, 27 (1980). 

^2 K-V Pham, M. Gabay, P. Lederer, Eur. Phys. Journ. B 9, 
573 (1999) 

G. V. Chester, L. Reatto, Phys. Lett. 22, 3 (1966); R. P. 
Feynman, Statistical Physics, (Addison- Wesley, New- York, 
1972) 

F. D. M. Haldane, Phys. Rev. Lett. 47, 1840 (1981) 
The zero mode basis (phonons and U± operators) forms a 
complete family of states: this was proved by Haldane in his 
computation of the grand canppical partition function of 
the Tomonaga-Luttinger modeB. Therefore Vq,j operators 
span the Fock space since Vq,j operators are products of 
Haldane's U± operators and of an exponential of phonons. 
Since the fractional states Vq_^ come from the decay of the 
Vq^j operators, they, along with the phonons span the Fock 
space. 

This is just the quantization of the superfluid phase: 
the boson operator is oc expjy'TFO and therefore 

^/n {0{L) — 0(0)) must be an integer multiple of tt. The 
total current around the ring is defined as the integral 
dx jnix) = Jr; therefore Jr = ^ (e(L) - 6(0)) and 

we may then write Jr — uKJ where J is an integer. 
''^ The previous equation assumes that both Q and n are pos- 
itive; if that is not so, it suffices to take their absolute value 
and replace the operators by their conjugate: for instance 



Q — ^ IQI and 



'^'^1/2 



w- 



1/2 



ix) 



D. Bohm, D. Pines, Phys. Rev. 85, 338 (1952) 
^® Recall that bosonization formulae for the electron opera- 
tors are obtained by making a Jordan- Wigner transforma- 
tion on the boson operators; likewise the fermionic (Q, J) 
selection rules {Q + J even) derive directly from the change 
in statistics: a (Q, J) excitation of the bosonic LL becomes 
a (Q, J' = J + Q) states due to the Jordan- Wigner trans- 
formation which attaches one fiux-tube to each particle {Q 
tubes to Q particles, which increases J by Q units) . Since J 
is even for the bosonic LL, Q -I- J' is even for the fermionic 
LL. 

S. C. Zhang, H. Hansson, S. Kivelson, Phys. Rev. Lett. 62, 
82 (1989); S. C. Zhang, Int. J. Mod. Phys. B 6, 25 (1992) 



Here and in the remaining of the paper for convenience 
we will consider the operators in real space and the wave- 
functions will be localized: the exact eigenfunctions are 
of course recovered by taking the Fourier transform of the 
localized (Wannier) states. 
*2 B. Horowitz, Phys. Rev. Lett 48, 1416 (1982) 

A. Luther, I. Peschel, Phys. Rev. B 12, 3908 (1975) 
The Jordan- Wigner transformation which maps spins onto 
fermions might lead one to think that the spectrum of 
the spin chain is identical to that of the Jordan- Wigner 
fermions. This is correct only for the spectrum of an even 
number of fermions added in (this corresponds to spin tran- 
sitions ASz ~ 2p where p is an integer): imposing periodic 
(PBC) or anti-periodic (ABC) boundary conditions forces 
the number of fermions to keep a given parity (even for 
PBC and odd for ABC), thus the mapping of spins onto 
fermions is not valid when one adds an odd number of 
fermions to the ground state. The spin fiip operator is not 
a Jordan- Wigner fermion since one also creates a fiux tube 
(the Jordan- Wigner string). This has led to some confu- 
sion in the literature. The ASz = continuum is correctly 
described by the excitations of the Jordan- Wigner fermion 
system since this corresponds to the addition of an even 
number of fermions. For the XX spin chain this means 
that the ASz = continuum is due to spin one states which 
are Jordan- Wigner fermions. But if the mapping were also 
valid for ASz = 1 excitations, this would mean that the 
continuum is also due to spin one states. This is not cor- 
rect because as we know this latter continuum is due to 
pairs of spinous. 

By adiabaticity, we mean a finite quasiparticles weight, as 
for the Landau quasiparticles which is an adiabatic defor- 
mation of the screened free electron. 
^'^ J. C. Talstra, S. P. Strong, Phys. Rev. B 56, 6094 (1997) 
I. Safi, PhD Thesis, Universite de Paris, Orsay, unpub- 
lished; ibid, Annales de Physique 22, 463 (1997) 
M. Stone, M. P A. Fisher, Int. Journ. Mod. Phys. B 8, 
2539 (1994) 

''^ K.-V. Pham, PhD thesis, Universite de Paris, Orsay, un- 
published. 

5" B. I. Halperin, Helv. Phys. Acta 56, 75 (1983) 
5^ The Tj matrix is distinguished from the K matrix to avoid 
a possible confusion because its eigenvalues are 1/Kc and 

It would be more correct to call this operator an anti-holon 
since its electrical charge is Qe = {—e)Q = — e; but apart 
from the sign of the charge all fractional particles consid- 
ered in this paper and their particle-hole conjugates have 
the same properties. So we call this state a holon. 
"'^ A naive decoupling of the electron operator would lead to: 

*<T = exp-iV7r(e<T - "I><t) 
= exp — tt/2 (Oc — <l?c) 
exp -i^/ -IT /2{Qs - $s) 

The first term is sometimes identified as the holon and it 
has a semionic statistics. The second term is likewise sup- 
posed to be the spinon. These identifications are in general 
incorrect because these states are unstable because of chi- 
ral separation (see the anisotropic spin chain). In the spe- 



26 



cial case of SU{2) symmetry, the above expression for the 
spinon is accidentally correct. 

C. S. Hellberg and E. J. Mele, Phys. Rev. Lett. 67, 2080 
(1991). 

M. Ogata, H. Shiba, Phys. Rev. B 41, 2326 (1990) 
F. Woynarovich, J. Phys. A 22, 2115 (1990) 
K. Penc and J. Solyom, Phys. Rev. B 47, 6273 (1993) 
The ambiguity which is rather technical is the following: 
there are four kinds of anomalous exponents (two chiralities 
and two modes). For a given excitation in a given mode, the 
sum of anomalous exponents for the two chiralities is fixed 
but their difference is not. The dimensions are therefore 
only known up to a constant; to fix their values, Prahm 
and Korepin make the unproven assumption that they are 
perfect squares if the excitation corresponds to a primary 
field. 

Note that in Frahm and Korepin' notations = Ni is not 
the spin density. 

The charge matrix and its inverse have the same properties, 
Z2 symmetry or SU (2) symmetry in the relevant cases. 
P. W. Anderson, Physics Today, 50, n° : 10, pp 42-47 
(1997) 

W. P. Su, J. R. Schrieffer, A. J. Heeger, Phys. Rev. Lett. 

42, 1698 (1979) 

P. Fendley, A. Ludwig, H. Saleur, Phys. Rev. Lett. 74, 3005 
(1994) 

M. Grayson et al, Phys.Rev.Lett. 80, 1062 (1998) 



27 



FIG. 1. Spectrum of the gaussian model with a bosonic 
Fock space in the charge Q sector. The energy at zero mo- 
mentum is A((5) = -KuQ^ / (2L) as a function of Q the number 
of particles added to the system. The spectrum for fermionic 
Fock spaces is identical if Q is an even integer. The contin- 
uum is enclosed within the straight lines which are supported 
by the parabolic enveloppe 7r«(fe/feF)^/(2L). 
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FIG. 2. Spectrum of the gaussian model for fermionic 
Fock spaces in the case of charge sectors for which Q is an 

odd integer. Notice that the energy has now local minima 
for momenta k = ±1,±3, ±5, ... (in units of kp) instead of 
fc = 0,±2,±4,... 
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FIG. 3. Spectrum of the XXZ spin chain for a trans- 
verse exchange J and anisotropy A; the continuum is en- 
closed within the curves E{k) = (7ra/2) | sin(fc) ] and 
E(k) = 7rasin(fc/2). The parameter a is related to the 
anisotropy A = cos (9 by a = sind/d. The linearized spec- 
trum found by bosonization is also shown. 
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